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Quantum error correction allows for faulty quantum systems to behave in an effectively error free 
manner. One important class of techniques for quantum error correction is the class of quantum 
subsystem codes, which are relevant both to active quantum error correcting schemes as well as 
to the design of self-correcting quantum memories. Previous approaches for investigating these 
codes have focused on applying theoretical analysis to look for interesting codes and to investigate 
their properties. In this paper we present an alternative approach that uses computational analysis 
to accomplish the same goals. Specifically, we present an algorithm that computes the optimal 
quantum subsystem code that can be implemented given an arbitrary set of measurement operators 
that are tensor products of Pauli operators. We then demonstrate the utility of this algorithm by 
performing a systematic investigation of the quantum subsystem codes that exist in the setting 
where the interactions are limited to 2-body interactions between neighbors on lattices derived from 
the convex uniform tilings of the plane. 



Quantum computers are a technological possibility be- 
cause there exist methods for building these comput- 
ers out of physical components that fail to operate in 
an error-free manner. The theory behind achieving this 
makes up the field of quantum error correction [IHS] and 
fault-tolerant quantum computing ["^JL^J . Of particular 
note is the threshold theorem for fault-tolerant quantum 
computing [HHTUl IT^ . This theorem says that if a quan- 
tum system decoheres slowly enough, and sufficiently 
precise control is maintained over the system, then ef- 
fectively arbitrary error-free quantum computations can 
be performed. The way that this is achieved is through 
the use of quantum information which is encoded across 
multiple quantum subsystems into a quantum error cor- 
recting code. 

Difi'erent quantum codes have different advantages and 
disadvantages for implementation in a fault-tolerant de- 
vice |13| . In this paper we undertake a study of an impor- 
tant class of quantum codes, quantum stabilizer subsys- 
tem codes [TlHTTj generated by measurements that are 
tensor products of Pauli operators. Part of the signifi- 
cance of this class of codes is that they can be used to im- 
plement passive fault tolerance by turning the measure- 
ment operators into interaction terms forming a Hamilto- 
nian that provide energetic protection against errors; the 
first example of such an approach was the toric code and 
related models due to Kitaev [T51 |TS] , and a plethora of 
related approaches have now been investigated [2m - [25| . 

Previous approaches for studying quantum subsystem 
codes have focused on using theoretical analysis to find 
and investigate new quantum subsystem codes. While 
powerful, theoretical analysis has some disadvantages: it 
is limited to the 'cleverness' of the analyst, and it can be 
prohibitively expensive to perform systematic searches 
of large parameter spaces to pick out the gems in the 
dust. In this paper, we present an alternative approach 
that uses computational analysis to accomplish the same 



goals. The advantage of this approach is that one be- 
comes limited by the power of the computer rather than 
the brain of analyslQ 

In this paper we develop an algorithm that computes 
the optimal subsystem code for a given set of measure- 
ments consisting of tensors products of Pauli operators. 
In the process of doing this we also develop a formalism 
that allows us to prove that the algorithm is correct and 
that the code it compute is indeed the optimal code for 
the given measurements. We also prove bounds on the 
running time of the algorithm that show that the algo- 
rithm terminates (relatively) quickly when the optimal 
code is not very robust to errors. Because of this prop- 
erty, the algorithm can be applied to sift through a class 
of possible measurements to determine which (if any) re- 
sult in a robust code. 

To demonstrate the use of this algorithm, we focus 
on classes of measurement operators where each mea- 
surement is limited in action to two qubits — that is, 
to operators taking the form Pi ■ Qj, where Pi and Qj 
are Pauli operators acting on respectively the i^^ and 
j'^ qubit of the system; examples of previous subsystem 
codes that have been constructed with this structure are 
the quantum compass model subsystem code [21] (in- 
cluding generahzations |29ll30]) and topological subsys- 
tem codes [31 . In particular we focus on systems where 
the measurement operators only couple qubits that are 
neighboring on a periodic lattice arising from the convex 
uniform tilings of the plane. We perform a systematic 
study of the codes on lattices arising from nine of the 
eleven such tilings, and present the results of this search. 



Of course, this is also the main disadvantage of this approach. 
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I. INTRODUCTION 

We begin by a brief review of the notion of quantum 
error correcting codes and in particular the subsystem 
stabilizer codes [l4| . 

In quantum computation we seek to reliably store and 
manipulate quantum information. Unfortunately, real 
quantum systems are open systems that couple to their 
environment and quickly lose their coherence through 
the process of decoherence. Even more troubling, when 



one wishes to manipulate quantum information one can 
only do this with a fixed precision. While considerable 
progress has been made in finding systems with long 
coherence times, inevitably current quantum computers 
will fail before they achieve anything close to the amount 
of computation needed, for example, to break a public 
key cryptosystem |32| . However it turns out that one can 
generally repair damage to quantum information as long 
as one knows the form that the damage took. Further- 
more one can build a 'trap' — that is to say, a quantum 
code — that tricks nature into giving up the information 
about what damage has occurred to the quantum system. 

The nature of codes is that they separate the space in 
which our computation lives from the space in which the 
physical information is stored; that is to say, although we 
design our quantum circuits to operate on some Hilbert 
space of qubits each of these qubits does not directly 
correspond to a physical qubit, but rather there is some 
isomorphism that relates the entire Hilbert space to 
the Hilbert space of physical qubits, To distinguish 
between these two Hilbert spaces, we call the Hilbert 
space of qubits in whose terms the computation is ex- 
pressed the computational space (or logical space), and 
the space of qubits which have physically been built the 
physical space. Merely building an isomorphism between 
these two spaces is not enough to allow us to correct er- 
rors. For one thing, we need to add extra qubits to the 
computational space that contain a record of the damage 
that we can read out; thus, we shall say that the full com- 
putational space is := ^ X where the qubits that 
live in ^ have the role of keeping a record of the errors 
that have been introduced by the environment, and the 
qubits that live in =S are the qubits in whose terms our 
quantum algorithm is expressed. 

We have to pick a strategy for reading out the infor- 
mation in ^ about the errors that have occurred on our 
system. One natural choice is to perform a single-qubit 
Pauli Z operator measurement on each qubit on M. In or- 
der to build the 'trap' element into our system, we need 
to ensure that whenever nature strikes at the physical 
space ^ and produces errors in a form that we intend 
to correct, this action must be isomorphic to a strike on 
the computational space that leaves a measurable record 
in For our choice of measuring Pauli Z errors, these 
are errors that are isomorphic to any operator that anti- 
commutes with a Z operator of at least one of the qubits 
in Note that although we speak of measuring the 
qubits in the measurement operator of interest in ^ 
is mapped to an operator in the physical space this 
isomorphic operator is referred to as a stabilizer, and the 
full set of operators on ^ which are isomorphic to our 
chosen measurement operators on ^ are referred to as 
the stabilizers of the code. 

Up to this point, the formalism we have described 
is known as stabilizer codes (6] l33H35j and its essential 
characteristic is that in determining the syndrome of the 
physical error, one makes a measurement of all of the 
qubits in What if, however, we relaxed this constraint 
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and only measured some of the qubits in Sil That is to 
say, what if we spht the qubits in into two categories: 
stabilizer qubits whose states we care about and which we 
measure to obtain an error syndrome, and gauge qubits 
whose states we do not care about. (The latter get their 
name from the fact that they provide a 'gauge' degree 
of freedom, i.e. a degree of freedom that is irrelevant to 
us.) Then we would have that ^ = =5^ x where ^ 
is the space in which the stabilizer qubits live, and is 
the space in which the so-called gauge qubits live; such 
a scheme is known as a stabilizer subsystem code |14| . In 
this case, we shall use the term stabilizers to denote the 
set of operators in ^ which are isomorphic to our chosen 
measurement operators of interest in y. 

At first there might not seem to be an advantage to 
this approach, since it essentially means adding qubits 
to our code that are 'wasted'; however, in practice sub- 
system codes have many advantages. The first advantage 
is that since we do not care about what happens to the 
gauge qubits, some quantum errors on the system will 
neither result in detectable errors nor destroy the infor- 
mation in the logical qubits [TlHTTl \7M A second 
advantage is that we no longer need our error-correcting 
measurements on the physical system to commute with 
each other, as long as they all commute with the sta- 
bilizers and logical qubit operators, since then the fact 
that they do not commute only affects the gauge qubits, 
which we do not care about [37]. This sometimes allows 
one to effectively measure a stabilizer which is a non- 
trivial fc-qubit measurement by using a series of two qubit 
measurements [3 7) . The individual measurements in this 
series do not commute (so they cannot be simultaneously 
measured) , however the stabilizer syndrome can nonethe- 
less be reconstructed from these measurements. A third 
advantage arises from the fact that subsystem codes of- 
ten require fewer measurements to diagnose errors than 
similar non-subsystem codes, which results in improved 
performance [131 I37j : counterintuitively, turning stabi- 
lizer codes into subsystem stabilizer codes often results 
in higher thresholds for fault-tolerant quantum comput- 
ing. Finally, subsystem codes can often be implemented 
in a more local manner than non-subsystem codes as ex- 
emplified by the quantum compass model code [HI [37] . 

There are now many examples of stabilizer subsystem 
codes in the literature. One of the first non-trivial sub- 
system codes to be described is a code related to the 
quantum compass model in two-dimensions [U] (33 [3^] ■ 
In the quantum model one considers a Hamiltonian on a 
two-dimensional square lattice where nearest horizontal 
neighbors couple the x component of their spins and near- 
est vertical neighbors couple the z component of their 
spin, so that the Hamiltonian is given by 

H = -A^{XijX^+i,j + ZjjZij+i), (1) 

where Pij represents the Pauli operator P acting on 
qubit at location (j,j). This model is interesting for a 
few reasons. The first is that the energy levels of this 



system can be best thought of as elements of a quan- 
tum error correcting subsystem code. The second reason 
is that the model provides some amount of protection 
from quantum errors because errors are energetically un- 
favorecj^ Many other examples of systems which have 
energy protecting properties are also known: the most 
famous being Kitaev's toric code in two and four-spatial 
dimensions 18, 19, 40j. The study of such systems is still 
in its infancy and one central question is whether there 
exist Hamiltonians with reasonable physical parameters 
(such as existing in three or fewer spatial dimensions and 
involving 2-body interactions [lUIll]) whose physics en- 
act quantum error correction on the system when the sys- 
tem is in contact with a thermal reservoir; such systems 
are called self- correcting quantum computers [241127) . In 
this paper we will talk about quantum subsystem codes 
from the perspective of active error correction where er- 
ror syndromes are identified through carefully engineered 
measurements, but it shall be understood that this for- 
malism can equivalently be seen from the perspective of 
passive error correction where errors are guarded against 
by carefully engineered interactions. That is, measure- 
ment operators in the active error correction picture are 
equivalent to interactions in the passive error correction 
picture. 

Because we ultimately want to build a system im- 
plementing our measurements, physical considerations 
typically constrain our measurements to be local, which 
means that they can be expressed in the physical space 
as a tensor product of single-qubit Pauli operators — i.e, 
for each measurement operator a we have that o := Pi 
where Pi is the pauli operator P acting on the i*'^ qubit. 
An important question then is which sets of local mea- 
surements give rise to useful quantum error correcting 
subsystem codes. 

Approaches to answering this question typically in- 
volve applying theoretical analysis with varying degrees 
of cleverness. In this paper we present an alternative 
approach. In section [ll| we present an algorithm which 
for every set of local measurement operators computes a 
quantum subsystem code that arises from the algebra of 
these operator^ Along the way we develop a formalism 
that allows us to prove not only that this algorithm is 
correct, but also that the code that it computes is opti- 
mal in the sense that there exists no other code arising 
from the same set of measurements for which the dis- 
tance of any of the logical qubits has been increased. 
This property makes this algorithm useful for analyzing 



^ Unfortunately, in this particular system the protection vanishes 
as the size of the lattice goes to infinity 1391 . but for small lattice 
sizes there is some protection from errors due to the energy level 
structure of the system 1381 . 

^ The code that we find is almost never unique, since among other 
transformations one can multiple every gauge and logical qubit 
operator by an element from the stablizers and end up with an 
equivalent code. 
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the properties of codes arising from measurements that 
are too overwhelming to analyze by hand. 

We shall also show that an important property of this 
algorithm is that it terminates (relatively) quickly when 
the distance of the code is small, which allows it to be 
used not only to solve for individual codes, but also to 
search through entire classes of sets of measurements to 
see if any have high-distance qubits. Motivated by pre- 
vious results demonstrating the utility of codes imple- 
mented using systems on a lattice, we undertake a sys- 
tematic investigation of codes where the measurement 
operators are restricted to the 2-body interactions aris- 
ing from the edges of periodic lattices derived from the 11 
regular tilings. In section [lll| we discuss our approach for 
applying the algorithm to perform a systematic search 
for codes that can be implemented on these tilings, and 
we the present numerical results that we obtained. In 
section [IV| we present our conclusions. 



A. Notation 

In this paper we adopt the following conventions for 
notation: 

• sets are denoted by a variable with a tilde, e.g. A; 

• sequences are denoted by a variable with an arrow, 
e.g. A- 

• operators and integers are denoted by using lower- 
case letters, e.g. a and i\ 

• collections of operators and pairs of operators are 
denoted by using upper-case letters with either a 
tilde or an arrow above them, e.g. O and O; 

• collections of integers are denoted by using lower- 
case letters with either a tilde or an arrow above 
them, e.g. k and fc; and 

• collections of other kinds of objects are typically 
denoted by capital letters in a fancy script. 



II. THEORY 

A. Construction of the subsystem code 

Remark This subsection describes by way of a construc- 
tive proof how to compute given a set of measurement 
operators the quantum code that can be implemented by 
these operators. For a listing of pseudo-code that imple- 
ments the algorithm described in this proof, see Table |l] 
near the end of this subsection. 

Although conceptually a subsystem code is an isomor- 
phism T such that ^«^^x^x,2 — that is, an 
isomorphism between the physical space of qubits and 
the computational space of qubits in whose terms our 



computation is actually expressed — we do not need to 
actually construct this isomorphism in order to be able 
to use the code. Since all of our work will be done on the 
physical system anyway, it suffices to know the operators 
in the physical space ,0^ that are isomorphic to the qubit 
measurement operators of interest in the computational 
space .5^ X ^ X =S, and it is exactly the operators on ^ 
that the algorithm we present shall computf]^ 

When one wants to define a qubit in terms of its mea- 
surement operators, it suffices to define two operators 
that anti-commute with each other but which commute 
with all of the others measurement operators that have 
been defined, since this gives us the X and Z measure- 
ments on the qubit which are sufficient to generate the 
full Pauli group (minus phases). Since working with such 
pairs of operators shall be a common theme in this algo- 
rithm, we shall introduce the following definition in order 
to simplify the language used to describe them. 

Definition A pair of operators is a conjugal pair in rela- 
tion to the set X when each of the operators in the pair 
commutes with every operator in X except for its conju- 
gal partner — that is, the other operator in the conjugal 
pair — should its conjugal partner be a member of X. 

Note that we have explicitly not required that the op- 
erators in the conjugal pair be members of X in order 
to be a conjugal pair in relation to it. However, should 
both operators be members of X, then neither operator 
can belong to a different conjugal pair with respect to 
X, since in that case there would be an operator in X 
(namely, its original conjugal partner) with which it anti- 
commutes that was not its conjugal partner in the new 
pair, leading to a contradiction. 

For convenience, we introduce the following additional 
definitions: 

Definition 1. *p is the group of Pauli operators — 
that is, the group of tensor products of the (unnor- 
malized) Pauli matrices — acting on the physical 
space , modulo phases] 

2. 'P(S') is the power set of 5, i.e. the set of all subsets 
of S ; and 

3. C&{S) is the centralizer of S*, that is the subgroup 
of elements in which commute with S] 

4. the function ^ : :P(*P) -> Vi^) is defined such that 
G{S) is the set of all possible products of operators 
in S" — that is, it is the set generated by S. 

We now introduce the main theorem of this subsection. 



* If one really wanted to, one could explicitly construct the isomor- 
phism ,"7 from these operators by computing the unitary operator 
which simultaneously diagonalizes a the maximal subset of com- 
muting measurements from this set of operators on ,53^, but in 
practice this is not particularly useful. 
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Theorem 1. Suppose we are given a sequence of Pauli 
operators, O. Then there exist sets of Pauli operators 
S C^, G C^, and L such that 

1. each of the operators in S U G U L is independent 
from the rest — i.e., no operator in this (unioned) 
set can he written as a product of other operators 
in the set; 

2. each operator in LU G is a member of a conjugal 
pair in relation to S U G U L; 



3. gCsijG)^g({d,})^ 

I andgiSUGUL) =Cq}(5) 



Remark This theorem follows, at least implicitly, from 
prior work on stabilizer codes [S] , the definitions of stabi- 
lizer subsystem codes given by Poulin [14], and the con- 
structive approach to finding such codes as exemplified 
in |24| . Because we wish to be constructive, however, we 
will present a full proof of this theorem and show how it 
gives rise to an algorithm for finding sets of Pauli opera- 
tors which satisfy Theoremjl] To be explicit, we note that 
S will be a set of stabilizers (or equivalently, generators 
for the stabilizer group), G will be a set of gauge qubit 
operators, and L will be a set of logical qubit operators 
(i.e., those on which the computation is performed). 

The main work in the proof of this theorem will be 
performed by proving several related propositions. First 
we shall show how the set G and a sequence S are con- 
structed from the sequence of operators O. Since we want 
our stabilizers to form an independent set of operators, 
we shall then show that through a Gaussian elimination 
procedure it is possible to extract a list of independent 
operators from a sequence S resulting in a set S. Finally, 
we shall show how using this same Gaussian elimination 
procedure we can transform a subset of the operators of 
S U G into a form that makes it trivial to compute the 
logical qubit operators L. 

Proposition 1. Suppose that we are given a sequence 
of Pauli operators O C *p. Then there exists a sequence 
of Pauli operators S C ^ and a set of Pauli operators 
G C <p such that 

1. all of the operators in S commute with each other 
and also all of the operators in G; 

2. each operat or in G is a member of a conjugal pair 
(Definition II A) in relation to {Si} U G ; and 



g [{s^} I 



uG =g({6,} 



Proof. Proof by induction. For the base case, note that 
if O is empty then 5^ := and G := trivially satisfy all 
properties. 

Now assume that the proposition holds for a sequence 
of length n — 1, and consider a sequence of operators O 
of length n. By the inductive hypothesis, we know that 
there is a sequence 5" and a set G' satisfying the proper- 
ties above for the subsequence of O consisting of the first 
n~l operators. Let a := On ■l\g^Q 0^ g^^oConiQ{g) — 

that is, the product of 0„ with the conjugal partner of 
every operator in G with which 0„ anti-commutes. This 
definition guarantees that o commutes with every oper- 
ator in G; furthermore, we can obtain 0„ back from o 
since every operator in G squares to the identity and thus 
On — o ■ Y[g(£G {O g}=o *''-'-'^jG('')' therefore we conclude 
that g ({S^} U G' U {o}) = g ({d,}) . 

If o commutes with every operator in S", then set 



S^ 



and G := G', and we are done. Otherwise, let s be some 

ferator in S" that anti-commutes with o, G :— G'U{s, o} 
SI' := /(S-), and S be the subsequence of S" with the 
identity operators removed, where 




s'-s {s',o} = 
s' otherwise 



Observe that by this definition, all of the operators in 
S commute with every operator in G, so property [l] is 
satisfied. Since the only difference between G' and G 
is the addition of s and o, which form a conjugal pair 
with respect to {Si} U G, we conclude that property pi is 
satisfied. Lastly, since s S G, we can form any operator 
in S' with products of operators in S and G, so therefore 
g U g) = ^ ({S^} U G' U {s, o}) = g ({Oj) , and 

so the final property is satisfied. 

We conclude by noting that since all of the operators 
in S and G were formed from products of operators in 
O, which are Pauli operators (i.e., members of the group 
*P), they are Pauli operators themselves. □ 

Remark A consequence of not requiring independence 



of 



3mark A consequence of not requiring independence 
the operators in O is that the operators S given by 



^ Here we use the notation {Oi} to refer to the set of elements in 
the sequence O. 



® Observe that neither o nor s can be present in G" since they 
commute with every operator in C?', so the new set G := G?' U 
{s, o} gives us a strictly larger set. This fact is irrelevant far 
as the proof is concerned, but it has the important consequence 
that a computer code implementing the algorithm described by 
this proof can append s and o to a list of gauge operators and 
assume that this list continues to form a set (i.e., a sequence 
without duplicates) without having to explicitly check for this. 
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Proposition [T] are not necessarily independent. Happily, 
since all of these operators can be expressed as tensor 
products of Pauli operators, we can construct a set of 
independent operators by performing an analog of Gaus- 
sian elimination. 

Proposition 2. Suppose that we have been given a 
sequence of Pauli operators which commute with each 
other, R. Then there exists 

1. a sequence S of n independent operators such that 

Q{{S^))=Q{{R^}), 

2. a sequence of n integers without duplicates in the 
inclusive range 1 . . .n, 



3. and a map p ; {1 . . . ri} — > {0, 1} such that Si 
only operator in S that anti-commutes with Pf 
where PP' :— Xk. and P|.^' :— Zk- 



the 



Proof. Proof by induction. For the base case, we observe 
that if R is empty, then the trivial sequences S :— ^ and 
k := % and the trivial function p : — > satisfy the 
requirements. 

Now suppose that we know the proposition holds for 
sequences of length — 1, and we are given a sequence S 
of length N . By our inductive hypothesis, we can apply 
the proposition to the first A'^ — 1 operators in R obtain 
sequences S' and k' of length n — ][^ and a map p' : 
— 1} — )■ {0,1} which all satisfy the respective 
properties of the theorem. Let 



R 



i=l...n-l 



n 

{««,P>f" }=0 



S'. 



We know that s commutes with every operator in S" 
because both s and every operator in S' are equal to 
products of operators in i?, which all commute with each 
other. Furthermore, since s is a product of Rjsf and a 
factor of S- for every i such that Rn and P|? anti- 
commute, and we know that S^ is the only operator in S" 
that anti-commutes with pj,? ^'^^ for i = 1 ... n — 1, it is 
therefore the case that s commutes with every member 
of the set {Pi.?'^'^'}j=i...„_i. Finally, smce s IS a prod- 
uct of -Rat and operators in S' , we can obtain Rn en- 
tirely from products of operators in {S-} U {s}, and so 

G{mu{s})^g(^m). 

If s is the identity operator, then let S :— S' and 
p :— p and we are done. Otherwise, we shall now show 



^ Note that n N in general, since some of the first A'^— 1 operators 
might not have been independent. 



that there must exist integers j € {!,..• ,N}\{k'^} and 
I £ {0, 1} such that s anti-commutes with pj'' , by demon- 
strating that if this were not the case then s would have 
to anti-commute with some element in 5', leading to a 
contradiction. 

Assume that s commutes with every operator in the 

set |p]'l ; j e {1, . . . , A^}\{fc^}, I e {0, 1}} . Recall- 
ing that s is a member of the Pauli group and thus 
a tensor product of single-particle Pauli spin matrices, 
and also that s commutes with every member of the set 

{pI^ '"^^^}i=i,,,n-i, we see therefore that s must be a prod- 
uct of elements from this set — that is, there is some 
subset ^ P C {P^y^'^^}i=i...n-i such that s = UoepO- 
However, from our inductive hypothesis we know that for 
every operator f G F there is an operator s' E S' that 
anti-commutes with / but commutes with the operators 
in P\{/}. Since s is therefore a product of a single op- 
erator that anti-commutes with s' and more operators 
that commute with s', we conclude that s and s' anti- 
commute, which contradicts our earlier conclusion that s 
commutes with every operator in S' . 

Now that we have shown that there exist integers j S 
{1, . . . , N}\{k'^} and I G {0, 1} such that s anti-commutes 
with p!'' , in terms of these integers we define 



S,. 



S- otherwise 



p{i) 



5' 

K l<i <n-l 
J i^n 

'p'{i) l<i<n-l 
I i — n 



and 



and we are done. 



□ 



Remark Proposition [2] is good for more than computing 
an independent set of generators from a commuting list 
of operators; it is also the key ingredient in computing 
the logical qubit operators. 

Proposition 3. Suppose that we have been given the ob- 
jects described in 1-3 of Proposition^^ Let S := {Si\i. 
Then there exists a set of operators L such that 

1. the operators in S U L are independent; 

2. every operator in L is a member of a conjugal pair 
with respect to 5 U i; 

3. g(^SUL^ = C^(^) — that is, the set generated 

by S U L is equal to the set of Pauli operators that 
commute with S. 
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Proof. Recalling that n is the number of elements 
in S (and S), let I be some sequential ordering of 
{1 . . . N}\{ki},, and then let L := {A;}, U {Bi}i where 

j = l...n 

To see that property [T] is satisfied, observe the fol- 
lowing. First, the operators in L are independent from 
the operators in 5* since none of them is the identity 
operator and they all commute with every operator in 
|pWO]|.^^ Second, they are independent from each 

other since for every i — 1 . . . |Z| we have that Ai is the 
only operator that anti-commutes with pjP^ and Bi is 

the only operator that anti-commutes with Pj^^^- Thus 

we conclude that all of the operators in S* U L are inde- 
pendent. 

Next, to see that property[2]holds, observe that for ev- 
ery choice of operators Ai and Sj we have (by intentional 
construction) that Sj either anti-commutes with two of 
the operators in the product forming Ai or none at all, 
and so [Si, Aj] = for alH = 1 . . . n and j = 1 . . . by 
the same reasoning we see also that [Si,Bj] = for all 
i = 1 . . .n and j = 1 . . .\l\. Furthermore, each operator 
Ai commutes with every operator in L except for its con- 
jugal partner Bi, since the only factor in Ai that could 
anti-commute with a factor contained within another op- 
erator in L is P;'.^' , and Bi is the only operator in L that 

contains a factor P^'*^' that anti-commutes with Xi. ; re- 
versing this argument, we also see that Bi commutes with 
every operator in L except for Ai. Thus, every operator 
in L is a member of a conjugal pair with respect to LUS. 

Finally, to see that property[3]holds, observe that since 
the operators in S commute they can therefore be si- 
multaneously diagonalized, which means that there is 
an automorphism on ^ that takes Si i~> p]^' for every 
i = 1 . . . n. The only operators that commute with ev- 
ery such pI^^ are those which do not contain any fac- 
tor of Pf' for i = l...n, and so C^j ({pf '},=!...„) = 

g {{PPh=l...n U {PP}^=n+l...N., i=0,l), which haS 2N - 

n generators. Since the automorphism preserves the 
number of generators in the centralizer, we thus conclude 
that C<^{S) has exactly 27V — n generators. Since S LI L 
contains independent operators which commute with ev- 
ery member of S, and furthermore |S'UL| ~ 2N — n, we 

thus conclude that g{s\JL^ = C^{S). □ 

With these building blocks in place, we now prove the 
main theorem: 



Proof of Theorem [7} By Proposition [Tj we know that 
there exists a list of operators S and a set of independent 
operators G satisfying the properties that are listed there. 
By Proposition!^ we know that there is an independent 
set of operators S that generate the same subgroup as S. 

Now let P be a maximal subset of commuting oper- 
ators in G — i.e., for each conjugal pair in G take one 
of the two operators — and then let O :— F U S. Since 
all of the operators in O commute, we apply Proposition 
[2] again to conclude the existence of the objects listed 
there, and then we immediately apply Proposition [3] to 
show that a set M exists with the properties listed there. 
We are not done yet, however, since there might be op- 
erators in G with which operators in M anti-commute, 
so we let 

L := {to • Y[ / ■ meM} 

{M,conjg(/)}=0 

where conj(=,(P) is the conjugal partner of F in the set G. 
This guarantees that the operators in L commute with 
every operator in S* U G, and so we are done. □ 

Remark A pseudo-code representation of the algorithm 
described by Theorem [l] is given in Table |l| 



B. Optimization of the logical qubits 

Remark A pseudo-code representation of the algoritm 
that will be described in this section is presented in Table 

Ell 

In general there are multiple sets of operators that sat- 
isfy the properties of [T] as is illustrated by the following 
Lemma: 

Lemma 1. Given conjugal pairs Q := (a, 6) and R := 

(c, d) in relation to some set X such that either a ^ c or 
h ^ d, we have that 

1. the pairs Q' :— (a • c,h) and R' :— {c,d ■ b) are 
conjugal pairs with respect to X\{Q, R} U {Q', P'}; 
and 

2. g {{a, b, c, d}) =gi{a- c, b,c,d- b}). 

Proof. 1. Since [a,c] — [a,d] — [b,c\ = [b,d] = 
{a, 6} = {c, d} — 0, we see therefore that [a ■ c,c] — 
[a-c,d-b] = [b,c] = [b,d-b] {a-c,b} = {c,d-b} = 0. 
Furthermore, since a, b, c and d commute with ev- 
ery operator in X\{Q, R}, so do a • c and d ■ b. 

2. Since b and c are Pauli operators and thus square to 
the identity, we have that a-c-c = a and d-b-b — d, 
and so g ({a, b, c, d}) — g {{a ■ c,b,c,d ■ b}). 

□ 
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1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 



Compute-Subsystem-Code(O) 

S^[] 
G^[] 
for o O 
do 



33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 



for {gx,gz) ^ G 
do 

if anti(o, gx) then o -f- 
if anti(o, gz) then o <- 

if o is identity then goto|3] 

for s ^ S 
do 

if anti(o, s) then goto 

goto|3] 

G ^ GU [(o,s)] 
i ^ 1 

for s' ^ § 
do 

if s' = s then goto 
if anti(s', o) 
then 

S[i] ^s' -s 

else 

S[{\ ^ s 
i ^ i + 1 

delete 



o- gz 
o- gx 
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16 



I ^ 
P 



28 call Gaussian-Elimination(S', 1, /, P) (Tahle\l]\ 

29 f^Syj[gx\{9x,gz)£G] 

30 call Gaussian-Elimination (f, |5| + 1, J, P) ("TafeZe^ 

31 L ^ 

32 for j <— 1 to number of physical qubits 
do 

if i G / then goto 

Ix^X, 
Iz ^ Z, 



32 



for {j,p,t) ^ iI,P,T) 
do 

if p = 
then 



if anti(t, Xj) then Ix ■'^ Ix ■ Zj 
if anti(f, Zj) then Iz ■'^ Iz ■ Zj 



else 



if anti(f, Xj) then Ix ■'^ Ix ■ Xj 
if anti(t, Zj) then Iz ■'^ Iz ■ Xj 

for {gx,gz) G G 
do 

if a.nti{lx , gz) then Ix ■'r- Ix ■ gx 
if a,nti{lz , gz) then Iz ■'^ Iz ■ gx 



Gaussian-Elimination(S, i, I, P) 

1 while i < |51 

2 do 



3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 



s ^ S[i\ 

for J to i — 1 
do 

{n,z)^{I[j\,P\j\) 
if z = 
then 

if anti(s, Xn) 

then s s ■ S\j] 

else 

if anti(s, Zn) 

then s <— s ■ S[j] 

if s is identity 
then 

delete S[i] 
goto[l] 

for n <— to number of physical qubits 
do 

if n G 7 then goto 

if anti(s, X„) 
then 

2 

goto 

if anti(s, Zn) 
then 



18 



29 



goto 



if 2 = 
then 



29 



else 



for j to i — 1 
do 

if anti(S[i],X„) 

then S\j] S[]\ ■ s 



for j to i — 1 
do 

if anti(S[i],Z„) 

then S[j] ^ S[j] ■ s 

append n to I 
append z to P 
S[i\ ^ s 

i -Ir- i + I 



50 return (5', G, L) 



TABLE II: Subroutine which performs a procedure 
analgous to Gaussian-ehmination on 5* to distih a set of 
independent operators from a possible dependent set of 
operators. This subroutine is cahed by the main 
subsystem code algorithm in Table |Tj 



TABLE I: Algorithm which computes the subsystem 
code generated by a given list of measurement operators 
O. The subroutine Gaussian-Elimination is listed in 
Table [n] 



As a result of this lemma, we see that we can take 
pairs of arbitrary conjugal pairs from sets G and L of 
Theorem [l] and replace them with different pairs per the 
recipe in Lemma [T] such that the properties of the theo- 
rem still hold. So given that these sets are not unique, 
the natural question is: What is the best choice of G and 
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L? To answer this, we observe that another criteria we 
would like for our code to satisfy is that it be as robust 
to errors as possible; in particular, we seek to maximize 
the difhculty of undetectable errors, which is defined as 
follows: 

Definition Given a set 5 C *p and operators I G *P and 
e e Ct^{S) which anti-commute (i.e., {l,e} = 0), we say 

that e is an undetectable error with respect to S acting 
on I. 

We assume that the 'difficulty' of an interaction between 
our physical system and its environment is related to the 
number of physical qubits in our system that are partic- 
ipating in the interaction. Thus, the natural metric for 
measuring the relative difficulty of an error is given by 
its weight, which recall is defined as follows: 

Definition Given an operator p £ *p — which recalls 
means that p is the tensor product of single-qubit Pauli 
unnormalizcd spin matrices — the weight of p is the num- 
ber of single-qubit operators in the product which are 
non-trivial (i.e., not the identity). So for example, the 
weight of / (g) / (g) / is 0, the weight of/0Z(8)/(8)Xis 
2, and the weight of Z (g) X (g) F is 3. 

For convenience, we introduce the following additional 
notation: 

Definition 

• the function 11 : P{^) — >■ *P is defined such that 
n(X) := n^ex ^ — that is, it is the product of 
the operators in X. 

• assuming we have a set of independent operators, 
Q, the function Gq : G{Q) 'P{Q) is defined 
(uniquely) such that for every q (z Q we have that 

9 = no6GQ(g) 

• the function w : *p — -jy is defined such that w{o) 
gives the weight of o; 

• the function Cg : (5) (Vo C<p) (S*) is de- 
fined such that eg{l) is the set of minimizers of w 
over the set : £ (^S^ , {o, = o| — that is, 

it gives the undetectable errors with respect to S 
acting on I that are of minimum weight; 



• the function M g is defined such that Mg | l is 



(P) is 



the function 



jV is defined such 



'S - 



w{o) for an arbitrarily chosen £ oe^ 



that I 

note that function is well-defined since all operators 
in the set have the same weight; 



the function m 5 : Crf, 



— ^ ^ is defined 



such that mg(Z, Z') := min{w^(^), aj^(Z')} — that is, 
it gives the smaller of the weights of the smallest 
weight errors acting on respectively L and L'; 



the sequence of |P| integers such that Alg IP) := 



(«); 



• the functions pi and p2 are defined such that, given 
(a, b) := X, we have that pi{x) :— a and P2ix) := b. 

• the function U : x ^) Vi^) is de- 
fined (for convenience) such that U := 

Uxex{Pi(^)^P'^(^)} — i^Sit is, it 'unpacks' a set 
of pairs of operators into a set of operators; in an 
abuse of notation, we shall also allow U to apply 

to sequences, so that U{P) := U (j^Pi}ij, and to 

individual pairs, so that if X is a single pair then 
U{X) := U{{X})- 

• finally, a choice of qubits stabilized by S, P, is a 
sequence of pairs of operators from the Pauli group 
such that 

1. no operator in U{P) appears in more than one 
pair in P; 

2. t/(P) CC^(^) 

3. every pair in P is a conjugal pair with respect 
to 5UC/(P); 

4. {uig opi)(Pj) = mg{Pi); and 

5. Mg{P) is an ordered sequence. 

Given the notation above, we now precisely define what 
we mean by the "best choice" of logical qubits. 

Definition An optimal choice of qubits stabilized by S is 
any choice of qubits, P stabilized by S, such that given 
any other choice of qubits, P', that is also stabilized by 
S and which satisfies {G o U){P) = {G o U){P'), we have 
that MiP), > M{P')i for aU 1 < i < |P| = |P'|0 

We now present an algorithm for computing the optimal 
choice of logical qubits from a set of input qubits. The 
key insight upon which the algorithm is built is that un- 
detectable errors acting on the space of logical qubits can 
never be eliminated entirely, so there will always be some 
operator on which they act. Thus, the goal of the opti- 
mization procedure is not to eliminate errors, but rather 
to contain them, so that they act on as few operators as 
possible. 

The optimization algorithm works by starting with an 
empty (and therefore automatically optimal) choice of 



* Note that since P and P' are sequences of conjugal pairs without 
duplicates they are therefore independent, and so if [Q oU)(P) = 
{Q o 0){P') then we know automatically that \P\ = \P'\. 
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qubits and a set of 'unoptimized' qubits, and making 
progress by gradually moving qubits from the unopti- 
mized set into the choice in such a way that preserves 
the optimality of the choice. The trick is that we want to 
delay as long as possible moving a qubit into the choice, 
until we have had every chance to improve it. Thus, we 
additionally keep track of a subset of pairs in the choice 
whose second members have yet to be used to contain 
an error, and then use them as much as possible to fix 
errors. That is, at every step in the algorithm, we scan 
for the minimal weight undetectable error acting on the 
set of operators consisting of both the second member of 
the pairs in this subset and all of the operators in the 
unoptimized set of qubits. If the minimal weight error 
acts on an operator in the first category, then we remove 
the pair from the subset and use this operator to fix this 
error wherever it occurs in both the second members of 
pairs in the subset and the unoptimized qubits. Other- 
wise, we pull out a qubit from the unoptimized set on 
which the error acts, use the first member in the pair to 
fix the error in the qubits remaining in the unoptimized 
subset, add the pair to the subset of qubits whose second 
members have yet to be used to contain an error, and 
then add it to the end of the choice. At this point our 
choice turns out to still be optimal because if there had 
been a way to make the qubit we just added any better 
by recombining it with other qubits in the choice then 
we would have already done so by now. 

This procedure is presented formally by means of the 
following inductive definition. 

Definition Let the function O be a map from a tuple of 
the form tuple {S, L) to a sequence of tuples each of the 
form {Q,P,s), where 

• 5 is a set of commuting Pauli operators; 

• i is a set of Pauli operators that are conjugal in 
relation to i U S"; 

• Q is a set of pairs of Pauli operators; 

• P is a sequence of pairs of Pauli operators; and 

• s is a sequence of integers from {0, 1} with the same 
length as P. 

The sequence is defined inductively. For convenience, 
we let the first index of this sequence be zero, and define 
0{S,L)o = (Z',0,0), where L' is the set of pairs such 
that U{L') and no operator appears in more than one 
pair in L', and is the empty sequence. Now assume 
that d{S,L)i is defined and that d{S,L), = {Q,P,s). 
If Q is the empty set, then 0{S, L)i is the last element 
of the sequence so that \0{S,L)\ = i + 1. Otherwise, 
d{S,L),+i := {Q',P',s'), where Q' , P' and s' are defined 
as follows. 

Let R := {p2{Pi) ■ i & {l...|P|},^ = 1} and X := 
U{Q) U R. Note that since Q ^% that therefore X 



Let 1^ be the minimal weight error with respect to S 
acting on any operator in X. There are two cases: either 
h acts on some operator in i?, or it doesn't and so must 
act on some operator in U{Q). The definition of Q\ P' 
and a! depends on which of these cases holds. 

Case 1: h acts on some operator in R 

Let k be the smallest index such that h acts on p2{Pk), 
and o '■=P2{Pk)- Define 



and 



7(a,6) 



X [h, x\=Q 
X ■ o {h, x} = 0, 



/ [h, a] = 
b {h,b}^0, 



Let s' Pi{Pk) ■ a ■ P where 
a n {{g{p2{P^).Pl{P^)) : i e {fc + 1 . . . |P|}, a, = 
and 

/3:=n({.g(a,6).g(6,a):(a,6)GO}). 

Then we define 

Q' ■.= {{f{a)J{b)):{a,b)eQ}, 

'Pi i < k 

(a', o) i = k 

P- ■■= <P^ i> k,di=0 

P, t> k,[h,p2{P^)]^Q 

. {pi{Pi),P2{Pi) ■ o) Otherwise 

_ jsi i^k 
1 i = k 



Case 2: h docs not act on some operator in R 

Let (J e Q be a pair such that h acts on one of its 
members, and without loss of generality assume that h 
acts on the first member since otherwise we can swap the 
members of the pair. Let o :— pi[q). Define 



/(^) 



X [h, x] = Q 
X ■ o {h, x} — 0, 



^ An observant reader may have noticed that we do not specify ex- 
actly how one goes about computing h. This was an intentional 
omission since the details are quite technical and fortunately they 
are irrelevent for proving that this algorithm works correctly as 
long as we can assume that h can be computed. Thus, the dis- 
cussion of how to compute h will be deferred until Section |ll B 3| 
when we analyze bounds on the running time of the algorithm. 



and 



I [h,a]=0 
b {h,b}^0, 



Let 



b" p2{q) ■ n ({g(a, b) ■ g{b, a) : (a, b) e Q\{q}}) , 



and 



6' 



b" [h, b"] = 
b" ■ o {h, b"} = 



Then we define 



{{f{a)Jib)):ia,b)eQ\{q}}, 

(P, I < \P\ 
\{o,b') i=\P\ + l 

\si i < \s\ 

1 1 t = \s\ + l 



Table [TTT| contains a listing of pseudo-code that uses the 
above algorithm to compute the optimal choice of qubits. 
For the sake of completeness, it includes additional steps 
that pertain to the details of how the minimal weight 
operator is computed, which will be discussed in more 
detail in Section Hi B 31 

In addition to proving that the above algorithm suc- 
cessfully constructs an optimal choice of qubits, we shall 
also provide a bound on its running time. In order to do 
this, we first need to precisely define what we mean by 
the running time for the purposes of this section. 

Definition We say that a computation can be performed 
in time x if the computation requires taking x products 
of Pauli operators. 

Of course, the number of products of Pauli operators 
is not the only metric that could serve as the gauge for 
the running time, but it suffices for our purposes. We 
now present the main result of this section. 

Tiieorem 2. Suppose we are given 

• a set of commuting Pauli operators, S , acting on 
N physical qubits; 

• a set of pairs, L C C(p(S'), conjugal with respect to 
U{Q) U 5; 

• and a set of Pauli operators C such that G{C) = 

then 0{S, L) is a sequence of finite length, and if 
{Q, P, s) is the last element in the sequence then P is 
an optimal choice of qubits such that Q{P) = Q{L), and 



Optimize-Logical-Qubits(5, G, L) 

1 iV^|L| 

2 fc ^ 

3 7n [true] * \ L\ 

4 nested function Query-Function(o) 

5 do 

6 for i to fc, {lx,lz) L[i] 

7 do 

8 if m[i] and anti(o, /z) 

9 tlien 

10 return (true, (1, i)) 

11 for k + lto N, {lx,lz) ^ L[{\ 

12 do 

13 if anti(o, /x) or anti(o, /z) 

14 then 

15 return (true, (2, i)) 

16 return (false, undefined) 

17 O ^ copy(S) 

18 for {gx,gz)^G 

19 do 

20 append gx and gz to O 

21 for {lx,lz)^L 

22 do 

23 append Ix and Iz to O 

24 P COMPUTE-PSEUDOGENERATORS(d) (Table VI) 

25 while k< N 

26 do 

27 (e, (c, I)) Find-Weight-Minimizer 



28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 

39 
40 
41 
42 
43 
44 
45 
46 

47 



(Query-Function, P) (Table 



{qx,qz) 
if c = 1 
then 



LI 



m[l] ^ FALSE 

for i <— Z + 1 to fc, (lx,lz) L[i] 
do 

if m[i] and {e, Iz} 
then 

qx <— qx ■ Ix 

L[i] <r- {lx,lz ■ qz) 

call FiX-LOGICAL-QUBITS 

{L,k,e,qz,qx) (Table 

L[l] ^ {qx,qz) 



IV) 



else 



if commute(e, gjf) 

then swap qx and qz 
{qx,qz) ^ L[l] 
L[l] ^ L[k] 
k^ k+1 

call Fix-Logical-Qubits 

{L,k,e,qx,qz) (Table\l 
L[k] {qx,qz) 



TABLE III: Algorithm which optimizes the logical 
qubits for a given subsystem code. 
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Flx-LOGlCAL-QuBlTS(L, k, e, qAyQB) 

1 for i ^ k to \L\, {lx,lz) ^ L[i] 

2 do 

3 if wciti{e,lz) and anti(e, /x) 

4 then 

5 qA QA ■ Ix ■ Iz 

6 L[i] <~ {Ix ■ qB,lz ■ qa) 

7 elseif anti(e,Zz) 

8 then 

9 qA qA ■ Ix 

10 L\i]^ {lx,lz -qB) 

11 elseif anti(e,Zx) 

12 then 

13 qA qA ■ Iz 

14 (/x ■ <3s,/z) 

15 if ANTl(e, gA) 



16 



then qA qA- qB 



17 return (qa, gs) 



TABLE IV: Algorithm which 'fixes' a subset of the 
logical qubits so that they are robust to a given error. 



furthermore it can he computed in a time that is in the 



set O \C\ 



10 



where d := M{P) 



.0 



The proof of this Theorem is rather technical and shall 
be split into several subsections. First we shall prove 
the existence of a condition that suffices to prove that a 
choice of logical qubits is optimal. Second we shall prove 
that the algorithm above constructs a choice satisfying 
this condition. Third we shall prove that the running 
time of the algorithm has the claimed bound. Finally 
we shall tie these results together to prove the Theorem 
above. 



1. Optimality condition 

How do we know that a choice of qubits is optimal? 
Intuitively, it should be sufficient to prove that a choice 
of qubits is optimal if we can show that there is no way 
that we can recombine qubits in the choice to form one or 
more qubits that are more robust than their component 
factors — that is, there is no way that any qubit can be 
"improved" by its involvement in such a product. This 
condition is stated formally in the following definition of 
unimprovable sets: 



A function / is said to be in the set 0{g) if / is asymptotically 
bounded by some fixed constant times g; formally / S 0{g) if 
and only if there exists constants c and xq such that f{x) < cg{x) 
for all X > xq. 



Definition An unimprovable set with respect to S is a set 
of Pauli operators, O, such that for any subset, X C O, 
we have that (wgo]^)(v^) — min^^ -^w^(x). We say that 

an unimprovable set O extends to Q if for all subsets 
X <Z duQ such that X n 6 ^ 
n)(X) < min^^j^^^ujgix). 



we have that {ujg o 



The following Theorem is the main result of this subsec- 
tion that proves that this condition is indeed sufficient to 
show that a choice of logical qubits is optimal. 

Theorem 3. If P is a choice of N logical qubits stabi- 
lized by S such that {pi{Pi)}i is an unimprovable set with 
respect to S that extends to U{P), then P is an optimal 
choice of qubits. 

Remark The intuition behind the proof of this Theorem 
is that because the set of first members of pairs is unim- 
provable and extends to the set of all members of pairs, 
we know that no qubit can be "improved" by recombin- 
ing it with one or more other qubits. Thus, the only 
way one could construct a better choice would by form- 
ing n + k independent qubits from products of n qubits 
(where k > 0), which intuitively should be impossible. 
Thus, we conclude that it is not possible for there to be 
a choice of qubits generated by the same qubits in this 
choice that is "better" than this choice. 

To assist us in proving this Theorem, we shall first 
prove a number of useful Lemmas and Propositions. We 
start with a simple Lemma that proves that taking a 
product of operators results in an operator that is no 
"worse" (with respect to its robustness to errors) than 
the worst operator in the product. 

Lemma 2. For any set of operators O, we have that 
(wg o n) (O) > mm^go (o) . 

Proof of Lemma. Any undetectable error with respect to 
S acting on (w^ o n)(<^) must also act at least one of the 
operators in O since otherwise it cannot anti-commute 
with the product. □ 

Remark In general, taking products of operators might 
result in an operator that is better than the worst oper- 
ator in the product because errors will cancel each other 
out — i.e., if two operators in the product anti-commute 
with an error then their product commutes with the er- 
ror. Thus, it is useful to state a condition under which 
we can be certain that this will not happen, so that the 
product is exactly as bad as the worst operator, which 
we do in the following Lemma. 

Lemma 3. Suppose we are given two operators a, 6 G 
such that LOg{a) < ujg{b); then ujg(a ■ b) = ujg{a). 

Proof of Lemma. Since i^g{ci) < ujg{b), there must be an 
undetectable error with respect to S that acts on a but 
not on b; thus, it must anti-commute with and hence act 
on the product a • 6, so that ujg{a ■ b) < ujg{a). Since 



11 Recall that M{P)^p^ is the distance of the best qubit in the ^^(^ . 5) > ^^(^^ LemmajS) we conclude that UJs.{a 



(optimized) code. 



b) 



,{a). 



□ 
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Remark Intuitively we should expect that it is not pos- 
sible to take n qubits and recombine them to form n + k 
independent qubits where fc > 0. To state this intu- 
ition in other terms, suppose we are given a set of con- 
jugal pairs C that are generated from some other set of 
conjugal pairs D. We know that every pair in C must 
have a member that includes a factor that is a first mem- 
ber of a pair in D (since otherwise the members of the 
pair cannot anti-commute), so let A be the set of first 
members of pairs in D. Our intuition then tells us that 
\C\ < \A\ = \I)\. The following Proposition states this 
fact formally: 

Proposition 4. Suppose we are given 

1. sets of independent Pauli operators Q and S; 

2. a non-empty set of conjugal pairs, C , with respect 
to U[C) U S, such that U{C) C g{Q); and 

3. a set A of independent Pauli operators with the 
property that for any conjugal pair X (a, 6) such 
that {a, 5} € Q{C), we must have that GQ{X)r\A ^ 



Then \C\ < \A\. 

Remark The basic idea behind the proof of this Propo- 
sition is that an analogue of Gaussian elimination can be 
used on the conjugal pairs to eliminate the presence of 
members of A from them; when we are done with this 
process, we can see that unless \C\ < \A\ we would have 
eliminated all members of A from some of the qubits, 
which contradicts the assumptions of this Proposition. 

The formal proof is somewhat technical and so we first 
introduce several Lemmas. First we prove a small helper 
Lemma that shows that it is possible to take a conjugal 
pair in which a given generator appears and rearrange it 
so that the generator only appears in the first member of 
the pair. 

Lemma 4. Let A :— {a,b) with {a,b} C Q{Q) be a con- 
jugal pair with respect to some set S, and a be some Pauli 
operator such that a G Gq{A). Then there exists a pair 
B := (c, d) such that 

1. {c,d}<zg{Q); 

2.oeGQ{c); 

4- (c, d) is a conjugal pair with respect to {^S\{a, 5}^ U 
{c, d}; and 

5. {gou){B) = igou){A). 



Proof. Let 



(c,d) : = 



{a,b) o£GQ{a),o(^GQ{b) 
(6, a) GQ(a),oG Gq(6) 

{a,b-a) o £ GQ{a),o £ GQ{b) 



Note that in any of the above cases, properties 1-3 are 
satisfied by construction, property 4 is satisfied because 
c and d are products of a and b which commute with 
every element in S\{a, b} and {c, d} = 0, and finally 
property 5 is satisfied because {a, b} C Q ({c, d}) and 
{c,d}cg{{a,b}). □ 

Remark This next Lemma contains the heart of this 
Proposition by introducing an analogue to a directed 
Gaussian elimination procedure. Specifically, it shows 
that if we have a generator a £ A that appears in one 
or more conjugal pairs, then we can take products of the 
conjugal pairs to eliminate it from appearing anywhere 
except in the first member of a single pair. 

Lemma 5. In the context of Proposition^ suppose we 
are given an element a £ A with the property that there 
exists a pair Y" £ C such that a £ G{Y"). Then there 
exists a conjugal pair Y and set of conjugal pairs D, all 
with respect to U (jY} U Dj U S, such that 

1. |^| = |C|-1 

2. {gou){{Y}\jb) = {gou){c); 

3. a£{GQop^){Y) butai{GQop2){Y); 
4- a,i {j^^f)GQ{D); and 



5. for every conjugal pair O 
GQ(O)ni\{a}^0. 



have that 



Proof. Proof by induction on the size of C . If C = {Y^"}, 
then apply Lemma |4] letting o := a, A := Y" , and Y := 
B, and we see that we have a pair Y which is conjugal 
with respect to {FjUS* and also such that {G°U){{Y}) = 
{goU){{Y"}). Let D := 0, and we see that the remaining 
properties hold trivially, so we are done. 

Now let us assume that this lemma has been proven 
for the case where \C\ = n — 1, and we are given a set C 
with n elements. Take any X" £ C\{Y"}, and apply the 
lemma to C\{X"}, A, a and Y" to obtain the objects 
Y' and D' described in this Lemma without the primes. 
If a ^ G{X"), then by the assumptions of the Lemma 
we know that G{X") n A\{a} ^ 0, so let Y Y' and 
D := D'U {X"}, and we are done. 

Otherwise, apply Lemma [4j setting A := X" , a :— a, 
and X' := B, and let X := (pi(X') ■ pi{Y'),p2iX')) and 
Y := {pi{Y'),p2{Y') ■P2{X')). Note X' and Y' are con- 
jugal pairs with respect to (j (^{X' ,Y'} U U S and 



14 



{X', Y'} n yDUSj =0, and so by Lemma |lj we con- 
clude that X and Y are conjugal pairs with respect to 
U (^{X,Y} UDjuS, and also that {G o U){{X,Y}) = 

[Q o U){{X' ,Y'}); since X' was obtained from apply- 
ing Lemma [4] to X" and a, we furthermore conclude 
that {g o U){{X,Y}) = {g o U){{X",Y'}). Since X' 
was obtained as a result of Lemma [4j we know that 
a e (G o pi)[X') but a ^ (Gg op2)(-'^'): we also 
know from the earlier recursive application of this Lemma 
that a e (Gopi)(y') but a ^ {Gq o p2){Y'). Thus, 

we observe that by construction, a £ [Gq °Pi){Y), and 

((GQop2)(r)UGQ(X)). 

Let D {X} U D' , and observe that \D\ = + 
1 = |G\{X"}| - 1 + 1 = |G| - 1, and also that {g o 

c/)({y} u ^) = (g o u){{x, Y} u = {g o c/)({x"} u 

({y'}U^')) = (^o^)({X"}U(G'\{X"}))^^ (^°^)(G). 
Furthermore, by the earlier recursive application of this 
Lemma we know that a ^ Gq{0) for every O E D' , 

so since we have also established that a ^ Gq{X), we 
conclude that a ^ Gq{0) for every O E D; since also 
know that every such O must also satisfy Gq(0) n A 7^ 0, 
we conclude that every such O satisfies GQ(0)n A\{a} 7^ 
0. □ 

Remark This next Lemma provides the small but im- 
portant result that we can always find a generator a that 
appears somewhere in the conjugal pairs; this has the 
consequence that we can now perform undirected Gaus- 
sian elimination (in contrast to the directed Gaussian 
elimination procedure described in the previous Lemma) 
by picking an arbitrary generator to eliminate rather 
than specifying a particular generator up- front. 

Lemma 6. In the context of Proposition^ there exists 
a Pauli operator a satisfying the assumption of Lemma 

m 

Proof. Take any pair Y' e G. By the assumptions of 
Proposition|4] we know that G{Y')nA ^ 0, which implies 
that there exists an element a € A such that either a e 
(Gq o pi){Y') or a G {Gq o P2)iY'). The existence of Y 

and D then follow immediately from the application of 
Lemma [5l □ 

Remark This final Lemma (inside the proof of Proposi- 
tion|4| shows using Gaussian elimination that there must 
be a number of generators from A present in the pairs in 
G that is equal to the size of G, since otherwise we could 
recombine the pairs in G to obtain a pair that includes 
no generator from A, contradicting the assumptions of 
Proposition |4j 

Lemma 7. In the context of Proposition^^ there exists 
a set of conjugal pairs, X , with respect to X \J S , and a 
subset of operators, O A, such that 



1. |X| = |0| = |G|; 

2. [g o U){X) = (g o U){C); and 

3. for every o € O , there is a conjugal pair Y E X 
such that o £ Gq{Y); 

Proof. Proof by induction. If G is empty, then the empty 
sets trivially satisfy this Lemma. 

Now suppose that we have proven this Lemma for 
|G| = — 1, and assume we have been given sets G 
and A such that |G| = N. Applying Lemma [g] to G 
and A we obtain the conjugal pair Y, the set of conjugal 
pairs D, and the element a described in the conclusions 
of that Lemma. Apply this Lemma recursively to the 
respective sets D and A\{a}, we obtain the sets X' and 
O' described (without the primes) in this Lemma; let 
X -.^ X'U {Y} and O := O' U {a}. Note that X is a set 
of conjugal pairs with respect to X U S since X' is a set 
of conjugal pairs with respect to X' U S, and we know 
that the operators in Y commute with every operator in 
every pair in X' since they commute with every operator 

in {g oU){D) ^ (g ou){x'). ^ 

First, observe that |d| = |0'| + 1 since a ^ 6' . Fur- 
thermore, a ^ {j^<^ij(x')G{x) since a ^ Uiec7(£i) G'Q(a:^) 
by Lemma g and {g o U){D) ^ {go U)iX') by recur- 
sive application of this Corollary. Thus, \X\ — \X'\ + 1 
since Y i X' as a e G{Y), and \X\ = |0| ^ \D\ + 1 ^ 

Second, observe that since {0 oU)(X') = (^o [/)(£)) = 
(g o U){C) by recursive application of this Lemma and 
{goU){{Y}UD) = (^ot/)(G) by Lcmmale) we conclude 

that {goij)(x) = {gou){{Y}ux') = {gou){{Y}uD) = 
{gou){c). 

Finally, observe that for every o e O we either have 
that o = a, in which case Y £ X and o e Gq{Y), or 

o G ^\{a}, in which case by recursive application of this 
Lemma we know that there is an operator Z £ X' Q X 
such that o e Gq(Z). □ 

Remark With the preceding Lemmas having performed 
the heavy lifting, the proof of Proposition [4] is quite sim- 
ple. 

Proof of Proposition Proof by contradiction. By 
Lemma [Tj there would have to exist a subset O C jl 
such that |G| = \0\ > \A\, which is impossible. □ 

Remark With the preceding Lemmas and Propositions, 
we now have all of the tools that we need to prove Theo- 
rem |3] Again, the idea behind this proof is that because 
the first members of pairs in the choice are contained 
in an unimprovable set, one cannot take products of the 
qubits in the choice in order to improve them; thus, the 
only way one could construct a better choice would by 
forming n + k independent qubits from products of n 
qubits (where fc > 0), which is disallowed by the result 
of Proposition |4] Hence, there can be no better choice. 
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Proof of Theorem^ Proof by contradiction. Let P' be 
some choice of qubits stabilized by S such that (G o 
U){P) ~ {G ° U)(P') (which automatically implies that 
1^1 = \P'\) ^'^d there exists some integer k such that 
Mg{P')k > Mg{P)k; in particular, let k be the smallest 
such integer, and let C := {P/ : i > k}. Let I be the 
smallest integer such that M§{P)i > Mg{P')k or |P| + 1 
if there is no such integer, and let A := {pi{Pi) : i>l}; 
note that since Mg{P')k > Mg{P)k we must have I > k, 
and hence \C\ > \A\. 

Take any conjugal pair O := {a,b) such that {a, 6} e 
G{C). Since a and b anti-commute, it must be the case 
that {pi{Pi)}i n Gjj^p^{0) ^ 0, because if every opera- 
tor in Gjji^p^ (O) were the second member of a pair in P 
then a and b would commute. Let c be a choice of a or 
b such that {pi{Pi)}i n Gq^p^{c) ^ 0. By Lemma jij we 

know that Mg{P')k < ujg{c) since c G G{C)- By the as- 
sumption of this Theorem that {pi(Pi)}i is an unimprov- 
able set that extends to U{P), we know that Mg{P')k < 
u}g{c) < min{ujg{x) : x € {pi{Pi)}i n Gjjf^p^{c)}. From 
these bounds we conclude that {pi{Pi)}ir]Gjjf^p-^{c) C A, 
and since {pi{Pi)}i n Gjj^p^{c) ^ we see therefore that 
^u{P)(^) n i ^ and so Gjj^p^{0) nA^dl. 

We have now demonstrated that for every pair O := 
(a, b) such that {a, b} e GiC), we must have G{0) nA^ 
0. Observe that this means that sets C and A match the 
descriptions in Proposition g (letting set Q {P,}i), 
and thus we see that it is impossible for |C| > \A\, and 
so we have a contradiction. We thus conclude that no 
such choice P' can exist. □ 



2. Correctness of the algorithm 

We now prove that this algorithm is correct — that is, 
that it terminates and outputs an optimal choice of log- 
ical qubits. We do so by proving the following theorem, 
which is the main result of this section. 

Theorem 4. Given a set of commuting Pauli operators 
S and a set of pairs L C C^(S') conjugal in relation to 

U{L) U S , the sequence 0(5, L) is finite and if {Q, P, s) 
is the last element then Q = and P is an optimal choice 
of logical qubits stabilized by S such that G{P) = G{L). 

Remark Before proving this Theorem, we shall first 
prove several related Lemmas and Propositions. 

Our ultimate goal is to expand the unimprovable set so 
that it includes at least the first member of every conjugal 
pair in the set of logical qubit operators, since this means 
that we have satisfied the optimality condition. Thus, 
we want to be able to add operators to this set while 
preserving the property of being an unimprovable set. 



The following Lemma shows that if we have an opera- 
tor o in a set X to which some unimprovable set extends, 
then if the smallest weight undetectable error acting on a 
acts on no other operator in X then we may move o to the 
unimprovable set to obtain a new unimprovable set that 
extends to X /{a}. The intuition here is that because said 
error acts only on o, it cannot be canceled by multiply- 
ing o by other operators, and so it is an "unimprovable" 
operator that can be included in our unimprovable set. 

Lemma 8. If O is an unimprovable set with respect to 
S that extends to X :— {o} U X' , and there exists an 
undetectable error, h, of weight ui g{o) that acts on a but 
not on any operator in X' , then O' := O U {o} is an 
unimprovable set with respect to S that extends to X' . 

Proof. Take any subset R C d'\JX' such that Rf^^' ^ 0. 
We need to show that uig (r) < min^g^p,Q,a;^(a) where 
r n(P). 

First consider the case where o ^ i?; in this case we 
have that R C OUX such that Pnd = PnO'^0, 
and so since O extends to X we conclude that (r) < 

™n„ei?,n6'^s(a) = mina6fln6''*^s(a} 

Now consider the case where PflO' = {o}. In this case, 
by the assumptions of this Lemma, we know h acts on o 
but not on any operator in X' which implies that h acts 
on r and so UJg{r) < w{h) — Wg(o) = jm-a.^^p^Q,uj g{a) . 

Finally we consider the remaining case where {o} C 
R n 6'. Let Z -.^ R\{o} ^ 0. Since 6 extends to 
X and PnO = Znd^0, we know that ijjg{x) < 
™n„ezno^s(y) ^'^ d < a;g(o), then w^(r) < 

d = mva^^pf^QiUj g(a) , and we are done. Otherwise, since 
d > LOg{o) we know that h acts on r since there can be 
no other operator in Z that anti-commutes with h, and 
so u:g{r) < w{h) < w^(o) = mm^^p^^,ujg{a) □ 

Remark In Case 1 of the algorithm we take an element 
that is a member of an unimprovable set and replace it 
with the product of this element times some elements in 
the set to which the unimprovable set extends. We want 
to show that this preserves the unimprovability of the 
set, and this is done in the following Lemma. 

Lemma 9. Suppose we are given an unimprovable set O 
with respect to S that extends to X . Let a be any element 
in O, and Z ^ O U X such that a ^ Z . Let a' := Tl{Z) 

and 6' (o\{o}) U {o'}. IfuJ-gi^d) ^ ujg{o), then O' is 

also an unimprovable set with respect to S that extends 
to X. 

Proof. Take any subset of elements i? C O'UX such that 
P n O' 7^ 0, and let x := U{R). We need to show that 
'^s(^) — ^"'^^^aeRnd'^si'^)- ^ then this follows im- 
mediately from the fact that RCOUX and Rf) O 7^ 
and O is an unimprovable set that extends to X, so as- 
sume that that 0' e R. Since 0' = I1{Z) and Z COUX, 
we conclude that the set T C O U X which is the sym- 
metric difference of Z and R satisfies the property that 
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X = n(T'). Note that o G T since o £ Z and o ^ O' and so 
o ^ R. Thus, TnO ^ 0, and so < min^gj,p|Qa;g(a) 

since O is an unimprovable set that extends to X. Thus, 
for us to show that ujg{x) < vam^^j^^Q,ujg{a), it suffices 
for us to show that vain^^f^^Qiijj g{a) = vain^^f^QUj g{a) . 

First note that since £^5(0') = oJ§{o), there is no el- 
ement z € Z O such that ujg{z) < (^§{0). Thus, 
any operator t E T O O such that u!g{t) < ojg(o) must 
also appear in i?n O', and vice versa; put another way, 
any operator that is less robust to errors than o must 
be present in both TnO and Rn O' together or nei- 
ther. Thus, if at least one such operator exists, then 
we conclude that min^g^y^^^, a;^(a) = min^^rp^QUJg{a) 
since in this case any minimizer of must be shared 
between the two sets. If no such operator exists, then 
since o E T O and o' G R O O' and there is no 
other operator present in either set with a smaller min- 
imum weight undetectable error, we conclude that o is 
the minimizer of over TnO and o' is the minimizer 
over R (1 O' and since ui§{o') = uig{o) we have that 

Thus we have shown that mhi^^j^^Q,ujg{a) = 

m.ui^^f^QUig{a), and since our choice of R was arbitrary 

we conclude that O' is an unimprovable set that extends 
to X. □ 

Remark In both cases of the algorithm we replace a set 

to which an unimprovable set extends with a new set that 
a product of elements in the old set. We want to show 
that the new set is also an extension of the unimprovable 
set, and this is proved by the following Lemma. 

Lemma 10. IfO is an unimprovable set with respect to S 
that extends to X, and X' is a set such that X' C Q{X), 
then O also extends to X' . 

Prool. Take any subset R! C O U X' such that A := 
R'dO ^ 0. We need to show that (a;) < min^g^a;^(t/), 

where x := n{R'). Note that since X' C g{X), there 

exists a set i? C X such that .t = 11 U , and so since 

O extends to X we conclude that Wg(a;) < min.y^^uig{y). 

Since our choice of R' was arbitrary, we conclude that 
that O extends to X'. □ 

Remark Most of the heavy lifting in this section is per- 
formed in the following Proposition, which uses induction 
to prove a number of properties about the output of the 
algorithm at every step. 

Proposition 5. Given a set of Pauli operators S and a 
set of pairs L C C^{S) conjugal in relation to U{L) U S, 

for every {Q, P,s) & 0{S, L) we have that 

1. g (U{Q) U U{P)) = g{L); 

2. U{Q) U U{P) C C^{S); 



3. 2(|Q| + |P|) = |Z|; 

4. \U{Q)i^lJ{P)\ = \L\; 

5. U{Q)f\U{P) = 0, and no operator appears in more 
than one pair in either Q or P; 

6. O is an unimprovable set of operators that extends 
to X; 

7. max^^oW§{o) < min^^^ix); 

8. {cOg opi)(q) = mg{q) for all q G P; 

9. M{P) is ordered; 

10. P is a choice of qubits stabilized by S; 

11. P is an optimal choice of qubits; 

where X := U{Q) U {p^iPi) ■ I < i < \P\,Si = 1} if 
P is non-empty and X := U{Q) otherwise, and O := 
U{P)\X. 

Proof. Proof by induction. It is easy to see that these 
properties hold for 0(5*, L)o = (L',0,0), so now as- 
sume that they hold for {Q,P,s) := 0{S,L)i, and let 
{Q',P',s') := 0{S,L)i+i. For convenience, define X := 
UiQ) U {p2{Pi) : 1 < i < |P|, Si = 1} (or X := U{Q) if P 
is empty), X' := U{Q') U {p2{Pl) ■l<i< \P'\,^i = 1}, 
O := U{P)\X and 6' := U{P')\X'. 

We now prove each of the conclusion above; note that 
in each conclusion we may assume that the conclusions 
prior to it have already been established, so we do so 
implicitly. 

Also, when we say that we are assuming we are in 
"Case 1" or "Case 2", we mean that we are assuming 
that {Q', P', s') followed from respectively Case 1 or Case 
2 in the definition of O. 

1. Examination of the definition reveals that Q' and 
P' are constructed entirely from products of ele- 
ments in g and P so that g (u{Q') U U{P')J C 

g{u{Q)uu{p)). 

If 0(5, L)j_|-i was defined using Case 1 then let 
(6, a) :— Pk and (6', a') := P'f. where k is the integer 
described in Case 1; otherwise let (a, b) q be the 
pair selected from Q in the definition and {a' ,b') 
be the last element of P' . Note that in either case, 
a = a'. 

In both cases, observe that for every operator € 
UiQ) U U{P)\{b} we have that either o or o • a 
is contained in U{Q') U U{P'), and since a is also 
contained in this set we see immediately that any 
operator in [/(Q) U [/(P)\{5} can be obtained from 
products of elements in U{Q') U U{P') (i.e., from 
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an operator in U{Q') U U{P') times possibly a). 
Thus we conclude that G (u{Q) U U{P)\{b}^ C 

Since b' is the product of 



g [U{Q') U C/(F') 

& with elements in G l^U (Q) U U {P)\{b} ^ and 
g(u{Q)UU{P)\{b}) C g(c7(Q')U(7(P')), we 
conclude that b e G {u{Q') U (P')) and so 

g (u{Q) u u[p)) c g (;7(Q') u ;7(p')^ 



Thus we have proven that g U{Q') U L/(P') 



g \U{Q) U ?7(P)j giL), and so we are done. 

2. This follows from the fact that L C 
C<p(^) ^ ^(L) C C^(^) and U{Q) U i7(P) C 

g (U{Q) U C/(P)) = giL) (which we just proved). 

3. By construction, either \Q'\ ^ \Q\ and |P'| = |P 
(in Case 1) or |Q'| = |Q| - 1 and |P'| = |P| + 1 (in 
Case 2). In either case we have that 2(|Q'| + |P'|) = 
2(|Q| + |P|) = |i|. 

4. Since the elements in L are members of conjugal 
pairs, they are therefore independent, and so we 
see that we need at least \L\ operators to gener- 
ate g{L). Thus we need \L\ < \IJ{Q') n [/(P')l < 
2(|(5'| + |P'|) = \L\, where the second inequality 
comes from the fact that a pair can unpack to at 
most two operators, and the last equality comes 
from the previous conclusion. We thus conclude 
that \U{Q')r\U{P')\ = \L\. 

5. By combining the previous two conclusions we see 
that \U{Q') U {7(P')I = 2(|Q'| + |P'|); if this con- 
clusion were false (i.e., an operator were repeated 
somewhere) then we would have that \U{Q') U 
U{P')\ < 2{\Q'\ + |P'|), which is contradicts our 
earlier results. 

6. First assume that we are in Case 1. Let k be the 
integer described in this case, (a, &) := P/j, and 
{a',b') := P'^. Note that a € O, and by con- 
struction a' := n(JL') where A' C O U X and 
A' f^O = {a}, and so since by the inductive hy- 
pothesis we know that O is a unimprovable set 
that extends to X we know that ujg{a') < w^(a). 
Since by the inductive hypothesis we also know 
that max^gQa;g(o) < min^gj^(x), by Lemma [2| we 
conclude that u!g{a') = u!g{a). Lemma ^ thus ap- 
plies to our situation and allows us to conclude that 

O" := (^0\{a}^ U {a'} is an unimprovable set that 

extends to X. Furthermore, si nce by construction 

allows us to con- 



there is an error of minimal weight that acts on b 
but not on any other operator in X', which means 
that by Lemma [sjwe conclude that 6" U {b} = 6' 
extends to X' . 

Now assume that we are in Case 2. Note that the 
only change from O to O' is the addition of a single 
element o that has an error that acts only on it 
but not on any other operator in X' . Note that 
since {o} U X' C g{X) we conclude from Lemma 
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X' C g{X) and b£X, Lemma 
elude that O" extends to {6}UX'. By construction. 



that O extends to {o} U X', and from Lemma [s] 
wc conclude that O' extends to X' . 

7. First observe that since X' C g{X), we con- 
clude from Lemma [2] that vtrni^^-^oj gix) < 

The difference between O and O' is the addition 
of a minimizer of ujg over X, o, and possible 
also the replacement of a single element. Since 

we conclude that since ujg{o) = min^gj^a;^(a) that 
therefore max^go^^^jW^(a) < min^g^,Wg(a). If 

O' = OU {o} then we are done. Otherwise, we are 
in Case 1 which means that we have also replaced 
an element in O; however, the operator we have re- 
placed it with is the product of an operator from 
O and operators from X, and since O is an unim- 
provable set that extends to X we conclude that the 
replacement can be no better than the operator it is 
replacing. Thus, max^gQ,aj^(a) < min^^j^, aj^(a). 

8. Since max^gQ,LLi^(a) < min^gj^.^,a;g(a), we immedi- 
ately conclude that (ujg opi)(P-') — mg{P[) when 

= 1. By the inductive hypothesis, we know that 
(a;^ opi)(P/) = m§{PI) where P/ = Pi] further- 
more, in both cases the pairs at the locations where 
di = are unchanged from P to P', and in each 
case this turns out to leave just a single location 
that wc still need to examine. 

In Case 1, this location is the index k described 
in that case, where — 1 and a'^ — 0. Since 
Pi(P^) is the product of a single element of O 
and elements from X, we conclude from the fact 
that O is an unimprovable set that extends to X 
that {ujg o pi)(P^) < {ujg o pi)(Pfc). By the in- 
ductive hypothesis we know that uiax^^QUjg{a) < 
mm ^^j^ujg (a). Because P2(Pfc) is a product of el- 
ements from X we conclude from Lemma [2] that 
^^^aex^si"-) ^ {^s°P-2)iPk)- Since pi(P) G O, 
we conclude that (w^ opi)(Pfe) < mm^^-^ujg{a). 
Combining all of these inequalities we reach the 
conclude that (aj^opi)(P^) < [ujg o p2){Pl) and 
hence (w^ o pi)(P(.) = mg(P^). 

In Case 2, this location is the end of the sequence 
P', but since the addition to the sequences is a pair 
of operators from X such that the first member is 
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a minimizer of ujg over X we conclude that (wj o 

9. By the inductive hypothesis we have that 
Mg{P)i = (uig o uig){Pi), and we have just shown 
that Ms(P')i = (wg o By the inductive 

hypothesis we know that Mg{P) is ordered, and so 
to prove that Mg{P') is ordered we need only check 

the places in the sequence where pi{Pi) ^ pi{P[). 
In both cases there is exactly one location where 
the first member of a pair is modified from P to 
P'. 

In Case 1, this is the index k defined in that 
which the first member was replaced with a prod- 
uct of that first member with elements in X. Since 
this member is in O, and since max^gQWg(a) < 
min^g j^a;g(a) (by the inductive hypothesis), we 
conclude from the fact that O is an unimprov- 
able set that extends to X that (w^ ° Pi){Pk) — 
{ujg opi)(Pfe), and so we conclude that Mg{P') is 
ordered. 

In Case 2, this is the end of the sequence P' 
where a pair was appended to P. Since the pair 
contains elements from X, and ma,x^^QUjg{a) < 

min^g j^a;g(a), we conclude that Alg{P')^p,^ > 
Mg{P)i for i < \P'\, and so we conclude that 
Mg{P') is ordered. 

10. The fact that P is a choice of logical qubits stabi- 
lized by S follows from directly from the previous 
conclusions. 

11. From the definition of an unimprovable set it is easy 
to see that since O' is an unimprovable set that 
extends to X' , it also extends to X'UO' ^ U{Q')U 
U{P'). Since {pi{q) : q e P'} Q 6' and U{P') C 
X' U O' , it is also easy to see from the definition 
that {pi{q) : q G P'} is an unimprovable set that 
extends to U{P') — that is, taking subsets does 
not affect the property of unimprovability. Thus, 
we conclude from Theorem [s] that P' is therefore 
an optimal choice of qubits. 

□ 

Remark Now that the heavy lifting has been done by 
the preceding Proposition, the proof of Theorem |4] is rel- 
atively simple. 

Proof of Theorem^ At every step in the algorithm, we 
either change an entry in s from 1 to or remove an 
element from Q. Since s is of finite length, as long as 
Q is non-empty there will be a step at which another 
element is removed from it. Thus, there is an index k 
such that if d{S, L)k = (Q, P, s) then Q = 0, and by 



definition this is the last element of the sequence. By 
Proposition [s] we know that P is optimal and also that 
0{P) = Q{L) (since Q is empty), and so we are done. □ 

3. Running time of the algorithm 

In this section we analyze the running time of the op- 
timization algorithm; the result is presented in the fol- 
lowing Theorem. 

Theorem 5. Suppose we are given 

• a set of commuting Pauli operators, S , acting on 
N physical qubits; 

• a set of pairs, L C C(p(S'), conjugal with respect to 
U{Q) U S; 

• and a set of Pauli operators C such that Q{C) = 

then the time needed to compute the sequence 0{S,L) is 
in the set O (^Cf + |Z|(|L| + d)3'' (^)) ' ^^^^^ 
A/(P)|^I and (Q, P, S) is the last element in the sequence 
(i.e., P is the desired optimal choice of qubits). 

Remark Before proving this Theorem, we shall first 
prove a number of related Lemmas and Propositions. 

The most complicated part of analyzing the running 
time of the optimization algorithm is analyzing the time 
needed to find the minimum weight undetectable error. 
In fact, the procedure for doing this was not even de- 
scribed explicitly in the algorithm, so we shall now ex- 
plain how we do it. 

The algorithm we employ is based on the Brouwer- 
Zimmermann search algorithm, which searches for the 
minimum weight binary string satisfying some property 
given a set of binary string generators endowed with a 
multiplication operation defined to be the exclusive-or 
operation. The Brouwer-Zimmermann algorithm works 
by using a Gaussian elimination analogue to express the 
generators in reduced row echelon form; it then performs 
its search by examining all products of r generators for 
increasing r. When all of the products of r generators 
have been enumerated, one knows that the set of strings 
that has yet to be enumerated has weight r-l- 1 or greater, 
since the row-echelon form means that every string has 
a column for which it is the only string with a 1 in that 
column, and so a product of k strings must have a weight 
of at least k. Thus, as the search proceeds, there is a 
growing lower-bound on the weight of the binary string, 
and the search halts when a string has been found that 
matches this bound. 

This algorithm cannot be immediately applied to the 
current problem because we are not working with binary 
strings, and in particular Pauli operators have a more 
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complicated multiplication operation that binary strings. 
Fortunately, in [43] White and Grassl showed that the 
Brouwer-Zimmermann enumeration can be generalized. 

The key difference between binary strings and Pauli 
operators is that binary strings only have two possible 
values in a given column, whereas Pauli operators have 
four. Thus, whereas we only need one element to gener- 
ate all of the possible values in a given column for a bi- 
nary string, we need two elements to generate all of the 
possible values in a given column for a Pauli operator. 
Thus, rather than working with generators, we instead 
work with a generalization that White and Grassl call 
pseudo-generators, which we shall define here as follows. 

Definition A pseudo-generator is a set of either 1 or 
2 Pauli operators. In an abuse of notation, we extend 
the functions U and Q to be respectively G ^ UgeG-9 
and G [Q o U){G) when applied to a set of pseudo- 
generators. 

Remark (The preceding definition does not follow that 
of White and Grassl exactly; it has been specialized to 
our situation for the sake of simplicity. ) 

Unlike 'normal' generators — i.e., Pauli operators — a 
product of generators is not a Pauli operator but rather 
a set of Pauli operators, which we define as follows. 

Definition Suppose we are given a set of r pseudo- 
generators G. Let X := {gig)\{I} : g € G}, Y he 
the r— ary Cartesian product of the r sets contained in 
X, and Z be the set consisting of the normal quantum 
operator product of the r operators in each r— tuple in 
Y. Then Z is defined to be the pseudo-product of the 
pseudo-generators in G. For convenience, we define a 
function 11 such that n(G) is the pseudo-product of the 
pseudo-generators in G. 

Remark The following Lemma places a bound on the 
size of the pseudo-product. 

Lemma 11. The pseudo-product of r pseudo-generators 
contains at most 3'' operators. 



IIB3 



has a 



Proof. Every set in X described in Definition . 
cardinality of either 1 or 3, and the size of Y is equal to 
the product of the sizes of all the sets in X\ since \X\ = r, 
we therefore conclude that the cardinality of Y and hence 
the number of operators in the pseudo-product is at most 
3''. □ 

Corollary 1 (to Lemma |11[) . Given a set of r pseudo- 
generators, G, then the set O := {/(o) : o € G{G)} can 
be computed in time 0({T + r)3^), where the time needed 
to compute f is in 0{T). 

Proof. From Lemma [TT] we know that there are at most 
3'' operators in the pseudo-product, so \0\ < 3^. Further- 
more, for every element in the set we first need to com- 
pute the corresponding operator in the pseudo-product. 



which requires r time since it is the product of r opera- 
tors, and then we need to compute /, which by assump- 
tion requires a time in 0{T). □ 

Remark In order to be able to place a lower bound 
on binary strings that have yet to be examined in the 
Brouwer-Zimmermann enumeration, we need the gener- 
ators of the binary strings over which we are searching to 
have the property that each generator has a column such 
that it is the only generator with a 1 in that column, so 
that products of r generators must have at least weight 
r. Because we want to similarly place a bound on unex- 
amined products of pseudo-generators, we generalize this 
property with the following definition. 

Definition A set of pseudo-generators G is said to be 
disjoint if for every g € G there exists some physical qubit 
k such that either Xk or Zk (or both) anti-conimutes 
with every operator in g{g)\{I}, but both Xk and Zk 
commute with every operator in Ug'eG\{g} Q{9')- 

Remark With the following Lemma, we show that the 
property of disjointness is exactly what we need to obtain 
the bounds that we want. 

Lemma 12. All of the operators in the pseudo-product of 
any r (distinct) pseudo- generators chosen from a disjoint 
set of pseudo- generators have weight of at least r. 

Proof. Every operator in the pseudo-product is the prod- 
uct of r factors, each of which is associated with some dis- 
tinct physical qubit k such that it anti-commutes with 
either Xk or Zk (or both) but every other factor com- 
mutes with both Xk and Zk] thus, the product must 
anti-commute with at least r single-qubit operators act- 
ing on distinct physical qubits, and so it must have a 
weight of at least r. □ 

Remark Now that we have the concept of a disjoint set 
of pseudo-generators and a result showing that an op- 
erator in a pseudo-product of r of them must have a 
weight of at least r, we present in the following Lemma 
an algorithm for searching through the space spanned by 
the pseudo-generators for an operator satisfying a given 
property. 

Lemma 13. Given a set of pseudo-generators G acting 
on N physical qubits and a test function / : *P — > {0, 1} 
such that /^^(l) n Q{G) 7^ 0, then a solution o such 
that f{o) — 1 and LOg{o) — min^i^g^^^^ /(o')=i^(''') '^'^'^ 

be computed in time O i {T + d)'i'^ (\^\ 



\w{o),\G\ 



where d 



and T is the time needed to compute f. 



Remark This Lemma follows directly from the results 
in (43|, though the proof is included here both for com- 
pleteness and also to show specifically how the results 
specialize to our case. A pseudo-code representation of 
this algorithm can be found in Table |V| 
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Proof. Define Cr to be the set of all operators such 
that if o S Cr then there is some subset of exactly r 
pseudo-generators from G such that o is contained in 
their pseudo-product. Note that UrCr — G{G), so for 
every operator o in the search space there is an integer r 
such that a E Cr- Corollary [l] shows that we can evaluate 
/ on every element of the pseudo-product of r pseudo- 
generators in time 0{{T + r)'S^), so since there are 

ways to choose r pseudo-generators from G we conclude 
that we can search Cr for a solution to /(o) = 1 in time 

O f(T + r)3'' 



From Lemma 12 wc conclude that w{o) > r for every 
o ^ Cr- By extension this means that w{o) > r for ev- 
ery o e Ut^lr ^r' , and therefore that if o ^ Ur'=^o ^r' 



Thus, if there exists an 
'«eU:rJoC',,,/(o)=i^(o) then we know 



and o £ Q{G) then w{o) > r. 
r such that r = min^ 

that r = min^ggj-Q^ /(o)=i''^("^) — that is, r is exactly 
the weight of the minimum weight solution to /(o) — 1, 
since any operator in the search space that isn't con- 
tained in Ur'l=^o ^r' must have a weight of at least r. Put 
another way, after having enumerated all of the elements 
in Ur'^^j check to see whether the smallest 

solution to / we have seen so far (if any) has weight less 
than or equal to r, and if so we are done since we have 
found the minimal weight solution. 

Now consider the procedure of searching through each 
Cr starting with r = 0. We know that we will eventually 
find at least one solution to /(o) = 1, since in this Propo- 
sition we have assumed that such an operator exists in the 
search space (by the assumption that f~-^{l)r\G{G) ^ 0). 
Furthermore, employing this procedure we will find the 
minimal solution a no later than after we have searched 
through Cr ior r — - - - w{o) , since at that point all of the 
unexamined operators have a weight greater than w{o). 
Thus, we conclude that we shall find the minimal weight 
solution after having searched at most all of the elements 

in Ur=o^'"'°'''''^'^ Gr, which we can do in time 



O 



min(-!«(o),|G|) 

E 

r=0 



iT + r)y 




C O ( (T + d)3'^ 



where d :— min(w(o), |G|). 



□ 



Remark The proceeding Lemma is rather general, so 
we shall show how it specializes to our case. First, 
however, we use the following three Lemmas to prove 
that our search space is generated by exactly N pseudo- 
generators. 

Lemma 14. Given a set of disjoint pseudo-generators, 
G, the largest subset X C U{G) such that X commutes 
has size \X\ < \G\. 
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Find-Weight-Minimizer(/, G) 

r ^ 1 
m oo 

while m > r and r < \G\ 
do 

for each H (- G such that \H\ — r, 
and each o in the pseudo-product of H 



do 



if weight (o) < m 
then 

(q^u) ^ f{o) 

if q is TRUE 
then 

m -f- weight (o) 
a <~ (o, It) 



if m = r 
then 



goto 



r + 1 
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return a 



TABLE V: Algorithm which finds the minimal weight 
operator in a given generating set that satisfies a given 
predicate. For the sake of convenience, we also allow the 
query function to return auxialiary information that is 
returned to the caller along with the minimal weight 
operator. 



Proof If X contained more that G operators then by 
the pigeon hole principle there would have to be at least 
two operators from the same pseudo-generator, and thus 
which did not commute □ 

Lemma 15. A set of disjoint pseudo-generators G acting 
on N qubits satisfies \G\ < N. 

Proof This follows directly from the definition and the 
pigeon hole principle. □ 

Lemma 16. For any set of commuting operators S act- 
ing on N physical qubits, if G is a set of disjoint pseudo- 
generators satisfying G{G) = C^{S) then \G\ = N - 

Proof- Since Ct^{S) contains a subset of N independent 
commuting operators, so must Q{G) and therefore tJ{G)] 
thus, Lemma 



14 



implies that \G\ > N, and combining 



this bound with that given by Lemma 15 we see that 
|G|=iV. □ 

Remark We now prove a Lemma which shows how the 
search specializes to the case of our qubit optimization 
algorithm. 

Lemma 17. Given 

• a set of Pauli operators S acting on N physical 
qubits, 
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• a set of disjoint pseudo- generator G such that 
g(G)^C{S), and 

• a non-empty set of Pauli operators Q such that QD 
5 = and every operator in Q is a member of a 
conjugal pair in relation to Q U S , 

then a minimal weight undetectable error acting 
on any operator in Q can be found in time 

O (^{\Q\ + d)3'^ (^^^^ where d:^ min{w{o),N). 

Proof. First observe that by Lemma [16] we know that 
|G|=iV. 

Define the function / : g{G) -> {0, 1} by 



1 3 q € Q such that {o, q} = 
otherwise, 



Note that solutions to / are undetectable errors acting 
on Q, and also that this function can be computed in 
time 0(|(5|) by checking the commutator for each ele- 
ment in Q. Furthermore note that for every operator in 
Q there is another operator in Q which anti-commutes 
with it, and also that Q C Ct]^{S). Thus, since Q is 

non-empty, there is at least one operator o € GiG) such 
that /(o) = 1. Thus, by Lemma 13 we know that we 
can compute a minimal weight solution to / in time 

O (^{\Q\ + d)3'' (^^^^ where d -.^ mm {w{o),N). □ 

Remark In order to make use of the preceding result, we 
need to have a set of disjoint pseudo-generators whose 
pseudo- pro duct covers our search space. However, we 
usually start instead with a set of ordinary Pauli opera- 
tors that generate this space. Thus, we shall now show 
that the former can be computed from the latter — i.e., 
that given a set of Pauli operators, we can compute a set 
of disjoint pseudo-generators that spans the same space. 
First we present a Lemma that provides a criteria suffi- 
cient to show that a set of pseudo-generators is distinct. 

Lemma 18. Given a set of pseudo-generators, G, if 
there exists a map p : U{G) — > {{Xk}k'^ {Zk\k) such 
that 

1. for every o € U{G), o is the unique operator in 
G{G) that anti-commutes with p{o) and 

2. for every g € G, the operators in g are both mapped 
by p to single- qubit operators acting on the same 
physical qubit k, and they are the only such oper- 
ators in Q{G) that are mapped by p to operators 
acting on k, 

then G is disjoint. 

Proof. For every g & G, we conclude from property 2 of 
p that there is some physical qubit k such that every op- 
erator in g is mapped by p to either or , and hence 



by property 1 this means that every operator in Q{g) 
anti-commutes with either or Z^. Since by property 
1 we know that the choice of or Z^. is different for 
each operator in Q{g), we conclude that if there is more 
than one operator in g then the product anti-commutes 
with both Xk or Zk ■ Finally, by property 2 we know that 
every operator in g' for g' € ~G\{g} commutes with X^ 
and Zk- □ 

Remark We now show that any set of operators that 
we are using to generate a search space can be expressed 
equivalently as a set of disjoint pseudo-generators. 

Lemma 19. Given any a set of Pauli operators, O, there 
exists a set G of psuedo- generators such that 

1. \G\<\d\, 

2. g(6) = g{G), 

3. the map p described in Lemma 
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exists for G, 



and G can be computed in time 0(|0p). 

Remark The structure of this proof bears some similar- 
ities to Proposition [2j In constrast with Proposition [2] 
however, in the setting of this Lemma we are working 
with operators that in general will not commute. 

Proof. Proof by induction. For the base case, we observe 
that if O is empty, then the trivial set G % and the 
trivial function p whose domain is the empty set satisfy 
the requirements. 

Now assume that this Lemma has been proven for sets 
of cardinality n — 1, and suppose we are given a (non- 
empty) set O of cardinality n. Take any operator o Cz O. 
By recursive application of this Lemma, we know that 
we can construct the set G' .— G and the function p' :— 
p described in this Lemma given O :~ 0\{o} in time 
0((n-l)2) =0(n2). 

Let 



n 



X. 



xec/(G'), 

{o,p(x)}=0 



Note that for every x G U{G'), it must be that o' com- 
mutes with p'{x), since o' is formed from a product 
that has either two factors that anti-commute with p' {x) 
(namely, o and p' {x)) or no operators that anti-commute 
w\i\ip'{x). If o' is the identity operator, then let G :— G' 
and p := p' and we are done. Otherwise, there must be 
some operator z g {{Xk}k U {Zk}k) /{p'{x) : x S U{G')} 
that anti-commutes with o. Define the function / by 



X ■ a' {x, z} — 0, 
X otherwise. 



and let G" := {{/(x) : x e g} : g e G'^ and p" := p' o 
f~^. Note that o' must be independent of the operators 
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in U{G'), because o' is not the identity and the product 
of o' with any subset of operators A C U{G') cannot be 
the identity since it must anti-commute with p'{a) for 
every a € A. Thus, / is a bijective map from U{G') to 
U{G") and hence is invertible, and so we conclude that 
p" is well-defined. Since, as previously discussed, o' com- 
mutes with p{y) for every y G U{G'), we conclude that 
multiplication by o' does not change whether any opera- 
tor X £ U{G') commutes or anti-commutes with p{y) for 
any y g U{G'), and so we conclude that the properties 
i that p' has in relation to G' (from the 
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listed in Lemma 
inductive hypothesis) are preserved in the transformation 
by / so that p" also has the same properties in relation 
to G" . Furthermore, since every operator x G U{G') was 
multiplied by a factor of o' if and only if it anti-commutes 
with z, we conclude that f{z) must commute with z, and 
thus every operator in U{G") must commute with z. 

There are two cases to consider: either there is no 
operator x e U{G") such that p"{x) acts on the same 
qubit as z, or there is exactly one, since if there were 
more than two then it would violate the properties of 
p" , and if there were exactly two then by construction 
o' would commute with z leading to a contradiction. In 
the first case, let G := G" U {{o'}}. In the second case, 

let G := {G"\{{x}}j U {{x,o'}}, where x is the single 

operator in U{G") such that p{x) acts on the same qubit 
as z. In either case, define 



p{x) 



p"{x) xeU{G"), 



observing that it is well-defined since U{G) = U{G") U 

W}. 

To prove conclusion 1, we note that G has at most 
one more element than G" and O always has one more 
element than 0\{o}, so conclusion 1 follows from this 
fact combined with the inductive hypothesis. 

To prove conclusion 2, we note that since o' (and 
hence o) is independent with respect to U{G'), then 
because of how G was constructed and the inductive 
hypothesis we have that g{G) = § (u{G") U {o'}) = 

g {u{G') u {o}) = g ((o\{o}) u {o}) = g{d)- 

To prove conclusion 3, we need to show that p satisfies 
the properties listed in Lemma [18] To prove the first 
property, we note that for every x G U{G) we have that 
either x G U{G"), in which case we have already shown 
that it is the unique operator that commutes with p{x) as 
this is true for the operators in U{G") as well as for o' (by 
construction), or a; = o', in which case this is still true 
since by construction o' is the only operator in U (G) that 
anti-commutes with z. To prove the second property, we 
note that due to the inductive hypothesis we need only 
consider the single change from G" to G, which consisted 
of either adding or replacing a pseudo-generator; in the 
first case (adding a generator), observe that we showed 



earlier that no operator y G G" is such that p"{y) acts 
on the same qubit as z, and in the second case (replacing 
a generator) , note that we added o' to the only generator 
in G" containing an operator y such that p"{y) acts on 
the same qubit as z; in either case, we see that the second 
property holds for G. 

Finally, we consider the running time. In addition to 
the 0{n?) time required to construct G', we required an 
additional 0(n) multiplication operations to construct o' 
and G"; hence the total running-time is O(n^). □ 

Corollary 2. Given any set of Pauli operators, O, there 
exists a disjoint set of pseudo-generators G such that 
|G| < \0\, g{0) — g{G) and G can he computed in time 

o(ion. 



Proof. This follows immediately from Lemmas 18 and 19 

□ 



Remark We now have the tools that we need to analyze 
the running time of the algorithm. 

Proof of Theorem First observe that from Lemma [18] 
we conclude that we can compute a set of disjoint pseudo- 
generators G such that g{G) — g{C) — Cf^{S) in 

time 0(|Gp). We will assume that this set of pseudo- 
generators is implicitly available to us throughout the 
algorithm so that we do not need to compute it more 
than once. 

At each step of the algorithm, we first need to find an 
operator, o, that has an undetectable error of minimal 
weight inside a set which we know from Proposition[5]has 



at most \L\ elements. By Lemma 
this operator can be found in time O 
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|L|+d)3'^ 



we conclude that 
N 
d 

where d :— mm{w{o), N) = w(o) (since the weight of 
any operator cannot be greater than N). After this 
has been found, examination of the algorithm reveals 
that the computation performed afterward takes a run- 
ning time in 0{\Q\ + \P\) = 0{\L\), where the equal- 
ity comes from Proposition [5] Thus, the total time 

needed for each step is in O y{\L\ + dji"^ 
O ( (\L\+d)3^'^ 



\L\ 



Let (Q' , P' , s*) be the second to last element of 0{S, L) 
and {Q, P, s) the last element. In the final step of the al- 
gorithm, we move the last remaining pair in Q' = {q} 
over to P', which means that the operator o with the 
minimal weight is a member of q. From Proposition [5] 
we know that that i^§{o) > mg{P[) for any i. Thus, 
at each step of the algorithm before this one we know 

that we spent a time in O ^(|L| + (i)3'^ (^^)) where 

d :— w{o) — i.e., a time no greater than the time 
spent on the last step. Since the algorithm requires 
at most \L\ steps we conclude that the total running 
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COMPUTE-PSEUDOGENERATORS(d) 

1 p^D 

2 i •(- 

3 while i < \0\ 



4 do 

5 o ^ d[i\ 

6 for J to i — 1 

7 do 

8 (n, 2) ^ 

9 if z = 

10 tlien 

11 if anti(o, X„) 

12 then o ^ o ■ d[j] 

13 else 

14 if anti(o, Z„) 

15 then o o ■ 0[j] 

16 if o is identity 

17 then 

18 delete d[i] 

19 goto|3] 

20 for n to number of physical qubits 

21 do 

22 ifanti(o,X„) 

23 then 

24 2 ^0 

25 goto [30 

26 elseif anti(o, Zk) 

27 then 

28 2 ^1 

29 goto [so 

30 if 2 = 

31 then 

32 for j to i - 1 

33 do 

34 if anti(d[fc],X„) 

35 then d[j] ^ d[j] ■ o 

36 else 

37 for j to i - 1 

38 do 

39 if anti(d[j],Z,0 

40 then d[j\ ^ d[j] ■ o 

41 append (n, z) to p 

42 d[i] •(- o 

43 i ^ i + 1 

44 G ^ 

45 for n to number of physical qubits 

46 do 

47 g^[] 

48 for i ^ to jdj ~ 1 

49 do 

50 (n',_)^p[i] 

51 if n = n' 

52 then append 0[i] to g 

53 ifff/D 

54 then append g to G 

55 return G 



TABLE VI: Algorithm which computes a set of disjoint 
pseudo-generators that generates the input set of 
operators. 



time, including that needed to compute the set of pseudo- 
generators, is in O (^\G\^ + \L\{\L\ + d)3'^ ( d^)) ' ^^'^^'^ 

ujg{o) = M(P)|P|, we conclude that d = M{P)^p^ (since 

there can be at most N qubits in the choice), and so we 
are done. □ 

III. PRACTICE 
A. Methodology 

In the previous section we presented an algorithm that 
computes the optimal subsystem code that can be im- 
plemented using a given set of measurements. The pro- 
cedure for optimizing the code requires an exponential 
amount of time, but fortunately the power of the expo- 
nential is a function of the distance of the best qubit in 
the code. Because of this property, this algorithm can 
be effectively applied to search over a set of choices of 
measurement operators to see if there is any good choice 
for implementing a code, since it can (relatively) quickly 
skip over the bad choices of measurements. 

In this section, we shall present an example of applying 
this algorithm to search for codes on quantum systems 
with the structure of a graph. That is, we assume that 
we have a system of qubits, 2-body Pauli measurement 
operators and a graph such that there is a bijection be- 
tween the qubits and vertices and between the edges and 
measurement operators, and also such that each mea- 
surement only acts on the two qubits corresponding to 
the vertices adjacent to its associated edge. Specifying a 
particular graph constrains the number of qubits and the 
types of measurement operators, but it still allows a great 
deal of freedom in the choice of the measurement opera- 
tor at each edge. In Figure [l] we illustrate an example of 
a graph with two possible such choices of measurement 
operator labelings; note that for the sake of generality we 
do not impose the constraint that the two operators in 
the 2-body measurement be identical. 

For reasons that will become clear, it turns out to be 
useful to specify choices of measurement operators in 
terms of ray labelings rather than edge labelings since 
the former is associated with vertices. Define a ray of a 
graph to be a pair consisting of a vertex and an edge ad- 
jacent to the vertex; note that every ray can be uniquely 
associated with an edge, and every edge can be associated 
uniquely with a ray for each of its two incident vertices. 
Thus, we can define a particular choice of measurement 
operators by labeling each ray in the graph with a single- 
qubit Pauli operator acting on the qubit of the incident 
vertex, and then letting the measurement operator as- 
sociated with each edge be equal to the product of the 
single-qubit operators in the edge's two rays. 

There is a natural symmetry of quantum codes that 
can be factored out to reduce the search space: the rel- 
evant properties of the code are invariant under single- 
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FIG. 1: A figure illustrating two possible labelings of a 

graph, which correspond to two possible choices of 
2-body measurements. In both graphs we see that there 
are four vertices and five edges, which indicates that our 

system is constrained to have four qubits and five 
2-body measurement operators. The edges (without the 
labels) indicate the pairs of qubits on which the 2-body 
measurement operators are constrained to act within 
our system. Within these constraints, we see in this 
figure two possible choices of measurement operators as 
specified by the two labelings of the edges of the graph. 



qubit rotations. That is, transformations such as swap- 
ping the X and Z operators at the location of a single 
physical qubit in every stabilizer, gauge qubit, and logi- 
cal qubit operator does not affect the code. Thus, when 
labeling the rays of a vertex, exactly which ray is labeled 
AT, y, and Z is not important; what matters is which rays 
commute and which rays anti-commute. We see therefore 
that we need only search over the possible ways to divide 
the rays into three indistinguishable groups, so that a 
vertex with n rays only has 1 + - — ^ — relevant labelings 
that need to be examined. 

The specific graphs we shall examine in this section 
are lattices generated by nine of the eleven convex vertex- 
uniform (also known as the "Archimedean" ) tilings of the 
plane — that is, those tilings with the property that ev- 
ery face is convex and every vertex has the same sequence 
of faces. Since these tilings have many translational sym- 
metries, we intentionally narrow our search to the set of 
labelings that share the translational symmetries of the 
latticcp^ Since the ray labelings must be preserved under 
these symmetries, we can partition the rays of the graph 
into equivalence classes such that two rays are equivalent 
if and only if they are related by a translation symmetry; 
thus we see that our narrowed search space is equivalent 
to the space of possible labelings of each class of rays in 
the lattice examined. Since there is a symmetry that can 
be factored out at each vertex (as discussed previously), 
we note that we can likewise partition the vertices into 
equivalence classes of vertices related by translation sym- 



Archimedean Tiling # Classes 

quadrille 1 

truncated quadrille 4 

snub quadrille 4 

isosnub quadrille 2 

hextille 2 

truncated hextille 6 

snub hextille 6 

dehille 1 

hexadeltille 3 

truncated hexadeltille 12 

rhombihexadeltille 6 



# 



Rays 
4 
3 
5 
5 
3 
3 
5 
6 
4 
3 
4 



# Labelings 
14 
625 
2,825,761 
1681 
25 
15,625 
4,750,104,241 
122 
2744 
244,140,625 
7,529,536 



TABLE VII: A table listing the number of vertex 
equivalence classes, the number of rays at each vertex, 
and the total number of labelings for each of the 11 
convex vertex-uniform tilings. By our scheme the 

/ 71-1 \™ 

number of labelings is equal to I 1 -I- - — ) i where 

m is the number of vertex equivalence classes and n is 
the number of rays at each vertex. 



metrics. If there are m vertex equivalence classes, and 
every vertex has n rays, then our search space consists 

( 



of 1 



total possible labelings. In Table 



VII 



^'^ This is not to claim that there are no interesting codes that 
break these translational symmetries; however, the investigation 
of such codes is outside the scope of this particular study. 



we list the eleven convex vextex-uniform tilings with the 
number of vertex equivalence classes, the number of rays 
at each vertex, and the total number of labelings. (Two 
of the eleven tilings, "truncated hexadeltille" and "snub 
hextille", had such a large number of possible labelings 
that we decided to exclude them from our search.) 

Note that we could furthermore refine our search to 
consist of those codes which also share the rotational 
symmetries of the lattice. We explicit avoid making this 
refinement because the existence of such codes as the 
quantum compass model code |24j indicates that there 
are good codes on lattices that require breaking the ro- 
tational symmetry of the lattice. However, we can use 
the rotational symmetries in a different way to reduce 
the search space as follows. Partition the labelings into 
equivalence classes such that two labelings are in the 
same class if and only if there is a rotational symmetry 
that relates them, and observe that all of the labelings 
in each class will give rise quantum codes with identical 
properties. Thus, we can reduce our search space to ig- 
nore redundant labelings by only examining one labeling 
in each equivalence class. 

Our search algorithm thus works in the following man- 
ner. We start by putting a total ordering on all of the 
lattice labelings (after having factored out the symmetry 
at each vertex.) We enumerate these labelings in order. 
For each labeling, we generate new labelings by applying 
each rotational symmetry to the current labeling. If any 
of these new labelings is less than the current labeling 
under our ordering, then we skip the current labeling be- 
cause we know that we have already previously examined 
an equivalent labeling. Although this algorithm proceeds 
serially through the search space, it can be parallelized 
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Archimedean Tiling 


# Non-redundant 


# Total 


Archimedean Tiling 


Symmetries 


Chosen 


quadrille 


10 


14 


quadrille 


p4m 


p4m 


truncated quadrille 


155 


625 


truncated quadrille 


p4m 


p4m 


snub quadrille 


706,881 


2,825,761 


snub quadrille 


p4g, p4, and pg 


p4 


isosnub quadrille 


743 


1681 


isosnub quadrille 


cmm 


cmm 


hextille 


11 


25 


hextille 


p6m 


p6m 


truncated hextille 


2392 


15,625 


truncated hextille 


p6m and p3ml 


p6m 


deltille 


58 


122 


snub hextille 


p6 


N/A 


hexadeltille 


594 


2744 


deltille 


p6m and p3ml 


p6m 


rhombihexadeltille 


904,741 


7,529,536 


hexadeltille 


p6m and p3ml 


p6m 








truncated hexadeltille p6m 


N/A 


TABLE VIII: A tabic listing the number of labelings 


rhombihexadeltille 


p6m 


p6m 



for each tiling that were not redundant under rotational 
symmetry transformations. This number was obtained 
by placing an ordering on the labelings and counting 
the number of labelings such that no symmetry 
transformation obtained a labeling less than the current 
labeling. Also listed for the sake of comparison are the 



total number of labelings from Table VII 



TABLE IX: A table listing the wallpaper symmetry 
groups (using crystallographic notation) for each of the 

11 convex vertex-uniform tilings, along with the 
particular group that we chose for our search. For two 
of the tilings no symmetry group was chosen because we 
decided not to search the tiling. 



by making use of n walkers, each of which starts at a 
different labeling (from to n — 1) and which proceed 
by examining the current labeling and then skipping di- 
rectly to the rfi^ labeling after the current one. In Table 
|VIII| we list the number of non-redundant labelings for 
each tiling. 

In order to preserve the rotational symmetries of the 
tiling, it is important that the lattice be constructed 
such that the center of the lattice is at a point of ro- 
tational symmetry. There is not a single unique center 
point that preserves all of the rotational symmetries of a 
given tiling, and furthermore for many tilings there are 
multiple rotational symmetry groups (known as "wallpa- 
per" symmetry groups) , each of which has a different set 
of center points. We thus chose the center of our lattice 
by picking the largest of the wallpaper groups present 
in the tiling and choosing the center to give rise to the 
rotational symmetries in that grou 



we 



In table |IX 

list the wallpaper symmetries for each of the 11 convex 
vertex-uniform tilings along with (where applicable) the 
particular symmetry group that we chose to utilize. 

As is usually the case in physical systems, it is impor- 
tant to pay careful attention to the boundary conditions 
of the lattice. In order to minimize boundary effects, we 
decided to put periodic boundary conditions on our lat- 
tices; care had to be taken to impose the periodic bound- 
ary conditions in such a way as to preserve the rotational 



Note that this approach does not mean that we have eliminated 
redundant labelings resulting from all of the symmetries in the 
lattice. For example, we have not eliminated labelings which are 
equivalent under rotations around a different point, nor which 
are equivalent under a glide-reflection symmetry. It is certainly 
possible to eliminate these labelings, but we choose not to do in 
this case in the interest of simplicity. An implication of this is 
that for many codes we expect to see many labelings giving rise 
to them that are equivalent under symmetry transformations but 
not eliminated by our approach. 



Tihng 


Maximum Radius ^ Qubits 


quadrille 


4 


64 


truncated quadrille 


6 


576 


snub quadrille 


5 


200 


isosnub quadrille 


8 


768 


hextille 


10 


600 


truncated hextille 


5 


600 


deltille 


8 


256 


hexadeltille 


3 


108 


rhombihexadeltille 


3 


162 



TABLE X: A table in which we show the maximum 
radius lattice that was completely scanned (i.e., such 
that every possible labeling was examined by the 
algorithm) for each tiling. To give a sense of the size of 
the lattices involved, we also list the number of physical 
qubits (corresponding to vertices) for the lattice with 
the maximum radius. 



symmetry group. For example, a boundary that only 
wraps from left to right and from top to bottom breaks 
some of the rotational symmetries for hexagonal tilings. 
In Figure [2] we illustrate how we placed the centers and 
the boundaries of the tilings. 

Due to limits on our computational resources, we were 
limited in the size of the lattices that we could search with 
the algorithm. We describe the size of the lattices using a 
quantity we call the 'radius', which is an integral quantity 
equal to the length of the lattice divided by the length of 
the smallest lattice defined for that tiling; the unit radius 
lattices are those illustrated in Figure figure:boundaries. 
In Table 1x1 we show the maximum radius lattice that was 
completely scanned (i.e., such that every possible labeling 
was examined by the algorithm) for each tiling. 

Since each labeling of every lattice results in a quantum 
code, we had to provide some criteria for our search al- 
gorithm to decide whether a code was interesting enough 
to log. We set our criteria relatively low: a code was 
deemed to be interesting if there was at least one logical 




(g) dcltillc (h) hcxadcltillc (i) rhombihcxadcltillc 

FIG. 2: A figure illustrating how we placed the centers and boundaries for each of the tilings that we scanned. The 
boundaries are periodic, so that edges that pass through one side of the boundary wrap around to the opposite side. 
Edges and vertices on a boundary arc merged with the corresponding edges and vertices on the opposite edge. Note 
that under this scheme there can be no vertices on a corner. For most tilings this will never happen, but it turns out 
that in the case of the deltille tiling there are vertices on the corners when the radius (smallest distance from the 
center to the boundary) is three times the radius of the unit cell; we thus ignore deltille tilings of these sizes. 
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qubit with distance three, that is if there was at least 
one logical qubit such that a single arbitrary error on 
that qubit can be corrected. This was done under the 
reasoning that as long as some of the logical qubits in 
a code are sufficiently useful to us to make implement- 
ing the code worthwhile, then we should not be troubled 
by the fact that there might be other logical qubits that 
are not useful because we can always ignore them (or, 
equivalently, classify them as gauge qubits). 



B. Results 

In the previous subsection we described the search 
space to which we applied the algorithm in order to 
computationally find possible codes that can be imple- 
mented using systems with 2-body interactions and a lat- 
tice structure following nine of the eleven convex vertex- 
uniform tilings. In this subsection we present the results 
of this search. The codes that we found are shown in the 
plots appearing in Figures |3] and |4] No plot appears for 
the deltille tiling because no codes were found for that 
tiling. It is worth emphasizing that these codes indicated 
in these figures are all of the (useful) codes that exist for 
the scanned lattices of that tiling given our constraints, 
since we scanned every possible labeling that was not 
redundant under a rotational symmetry transformation 
about the center. 

Observe that two kinds of trends appear frequently 
in the results: codes that grow in distance but remain 
constant in the number of logical qubits as the radius 
increase, and codes that remain constant in distance 
but grow in the number of logical qubits as the radius 
increases. The former trend appears in the quadrille, 
snub quadrille, isosnub quadrille, hexadeltille, and rhom- 
bihexadeltille tilingj^ The latter trend appears in the 
truncated quadrille, snub quadrille, hextille, truncated 
hextille, hexadeltille, and rhombihexadeltille tilings. In 
many of the tilings there are also codes that were found 
that do not seem to belong to an obvious trend. 

In the follow subsections we will focus on some specifics 
of the results for each of the tilings. 



1. quadrille 

For the quadrille lattice, we only saw one labeling, il- 
lustrated in FigurejS] that resulted in an interesting code. 
This labeling corresponds to the compass model code, 
and the algorithm correctly found that the distance of 
the code grows linearly with the radius of the lattice and 
is exactly equal to the square root of the number of qubits 



in the lattice. This result is not terribly surprising, but 
it is good to see that our search technique employing the 
algorithm can correctly duplicate known results. 



2. truncated quadrille 

There are three kinds of codes that appear in this tiling 
where the number of qubits increases with the radius: 
two where the distance is fixed at 4, and one where the 
distance is fixed at 3. For the best two of these three 
kinds of codes, the number of logical qubits (l) is related 
to the radius (r) by / = (2r — 1)^. Since the number of 
physical qubits (n) is given by n = (4r)^, the number 
of logical qubits per physical qubit is thus given by ^ = 

(I — ^)^, a quantity which converges to ^ as r — >■ co. 
There were four labelings with this property that we saw 
in our search: two with distance 3 qubits (illustrated in 



Figure 
Figure 



and two with distance 4 qubits (illustrated in 



3. snub quadrille 

There are two kinds of interesting codes found in this 
filing. First, we saw exactly one labeling that has the 
property that the distance is four and the number of 
logical qubits (l) is given by I = 2r{r — 1), where r is 
the radius of the lattice. Since the number of physi- 
cal qubits (n) is given hy n — 8r^, this means that the 
number of logical qubits per physical qubit is given by 



1 



I ^ 2r(r-l) ^ 1 
n Sr^ 4 

is illustrated in Figure [7] 



as r — >■ oo. This labeling 



Second, more usefully, we saw twelve labelings which 
result in a code that has one qubit whose distance grows 
with the size of the lattice. Two of these labelings are 
illustrated in Figure [Sj 



4- isosnub quadrille 

We only saw two labelings of the isosnub lattice that 
result in useful codes, both of which only have a single 
qubit that seems (assuming that the trend seen in Figure 
|4] can be extrapolated) to have a distance that grows 
with the radius of the lattice. These two labeling are 
illustrated in Figure |9j 



Note that where the former trend was present, the maximum 
radius that we scanned was often quite limited; this is due to the 
exponential explosion in the cost of finding the optimal code as 
a function of the distance of the code. 



5. deltille 

We scanned this tiling up to a radius of eight; no in- 
teresting codes were found in any of the 122 labelings. 
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Codes found for the 
quadrille tiling 




3 

(36) 



4 

(64) 



1 2 

(4) (16) 

Radius 

(Number of physical qubits) 



Codes found for the isosnub quadrille tiling 




1 2 3 4 5 6 7 8 

(12) (48) (108) (192) (300) (432) (588) (768) 

Radius 

(Number of physical qubits) 



(a) 



(b) 



Codes found for the truncated quadrille tiling 



(T 
O 

'So 
o 




1 2 3 4 5 6 

(16) (64) (144) (256) (400) (576) 

Radius 

(Number of physical qubits) 

(c) 



o 
'5b 

o 



(U 



3 

z 



Codes found for the 
snub quadrille tiling 




1 2 3 4 5 

(8) (32) (72) (128) (200) 

Radius 

(Number of physical qubits) 
(d) 



FIG. 3: A figure containing plots of the results from scanning the (a) quadrille, (b) isosnub quadrille, (c) truncated 

quadrille and (d) snub quadrille tilings. Every polygon in the plot corresponds to a code that was found with a 
distance equal to the number of sides of the polygon, so that triangles indicate distance three codes, squares indicate 
distance four codes, etc. The position along the x-axis indicates the radius of the lattice where the code was found, 
where the radius is an integer defined to be the length of the lattice divided by the length of the smallest possible 
periodic lattice for the tiling; it also indicates the number of physical qubits in the lattice where the code was found, 
which appears on the a;-axis just under the value of the radius. The position along the y-axis indicates the number 
of logical qubits with that distance in the code. Note that in cases where multiple codes were found for the same 
radius and with the same number of logical qubits, multiple polygons are drawn, so that for example in plot (c) we 
see several cases in which a distance 3 code (triangle) and a distance 4 code (square) were found that were in lattices 
with the same radius and also had the same number of logical qubits. 



hextille 



In this tiling we saw four labelings which resulted in 
two kinds of interesting codes: two of the labelings (illus- 



trated in Figure 10a I resulted in codes of distance 3 that 



were present for every value of the radius, and two of the 
labelings (illustrated in Figure 10b I resulted in codes of 
distance 4 that were only present for even values of the 
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Codes found for the hextille tiling 
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(24) 


1 

3 

(54) 


1 

4 

(96) 


— 1 — 
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1 

6 
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1 

7 

(294) 


1 1 
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(384) (486) 
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10 
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(Number of physical qubits) 
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Codes found for the 
rhombihexadeltille tilins 



'5b 
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3 
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16 
14 
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# 




13 
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1 
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1 
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(18) 


(72) 


(162) 
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Codes found for the 
truncated hextille tiling 




3 

CT 

'5b 
o 



S-H 



3 

z 



1 2 3 4 5 

(24) (96) (216) (384) (600) 

Radius 

(Number of physical qubits) 
(c) 



Codes found for the 
hexadeltille tiling 




1 2 3 

(12) (48) (108) 

Radius 

(Number of physical qubits) 
(d) 



FIG. 4: A figure containing plots of the results from scanning the (a) hextille, (b) rhombihexadeltille, (c) truncated 
hextille and (d) hexadeltille tilings. See the caption of Figure [s] for an explanation of how to interpret these plots. 
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x j xx i xx i xx l x 



x j xx i xx i xx l x 



x j xx i xx i xx l x 



x l xx i xx i xx l x 



FIG. 5: A figure illustrating the single labeling of the 
quadrille tiling (on a radius 2 lattice) that results in a 
useful quantum code. 



radius. The former resulted in codes which had a number 
of logical qubits (Z3) given by I3 — ilnlKlii^^ where r is 



the radius, and the latter resulted in codes which had 
a number of logical qubits {I4) given by I4 = r{r + 3). 
Since the number of qubits (n) is given by n = 6r^, 
we have that the number of logical qubits per physi- 
cal qubit for the distance and distance 4 codes were 
given respectively by da = ^ = ^^^d 
(i4 = ^ = i (1 + ^) ; as r — 7> 00, we have that ^3 — > 
and ^4 — > ^. 

It is interesting to observe that there is no dis- 
tance/qubit count trade-off in this tiling. As long as the 
radius is even, the distance 4 code is superior in both dis- 
tance and logical qubit count over the distance 3 code. 



7. truncated hextiUe 

In the truncated hextille there are four kinds of codes 
where the number of qubit increases with the radius: 
three with the distance fixed at 3, and one with the dis- 
tance fixed at 4. 

The best of the distance 3 codes has the number of 
logical qubits (Z3) given by I3 = 2r^ — 1, where r is the 
radius of the code. Since the number of physical qubits 
(n) is given by n = 24r^, the number of logical qubits 

per physical qubit is thus given ~ T2 ~ ^) 

quantity which converges to as r — > 00. The two la- 
belings we saw which give rise to this code are illustrated 
in Figure [TT] 

The best of the distance 4 codes has the number of 
logical qubits {I4) given by I4 — 2r(r — 1), and the number 
of logical qubits per physical qubit is thus given by — = 



(1 — i), a quantity which converges to as r — >■ cx). 
We see from this analysis that although the best distance 
4 code contains fewer logical qubits than the best distance 

3 code, they both converge to the same number of logical 
qubits per physical qubit in the large radius limit. One of 
the nine labelings we saw which give rise to this distance 

4 code are illustrated in Figure [12) 



8. hexadeltille 

There are many codes that appear in the hexadeltille 
tiling, but it is difficult to draw conclusions about trends 
due to the limit on the size of the lattices that were 
scanned. The good news, though, is that the reason 
why scanning larger radii was difficult is because there 
is a code in this tiling with a qubit whose distance grows 
with the radius of the lattice. One of the nine labelings 
that we saw with this property is illustrated in Figure 13 



9. rhombthexadeltille 

This tiling is interesting because it had many more la- 
belings that resulted in codes than all of the other tilings 
combined; specifically, for the rhombihexadeltille tiling 
we saw 48,807 labelings that resulted in useful codes, 
whereas for all of the other tilings combined we saw only 
421 labelings that resulted in useful codes. This is even 
more remarkable considering that the largest lattice we 
were able to scan for the rhombihexadeltille tiling was 
smaller than that the for most of the other tilings. 

As can be seen in Figure |4| this tiling is also interest- 
ing because it features so many different kinds of codes, 
including both codes that seem to grow in the number of 
logical qubits with radius and codes that grow in distance 
with size. It is the only tiling that features a lattice that 
contains a labeling resulting in a code for every distance 
up to 6. 

The rhombihexadeltille tiling is the only tiling we have 
seen which has code both with a distance greater than 4 
and multiple qubits; we saw six labelings which resulted 
in codes with distance 6 and two qubits, and four la- 
belings which resulted in codes with distance 5 and four 
qubits. In Figure[T4]we show an example of each of these 
labelings. 

Two of labelings resulted in codes with the highest 
number of qubits - 16 logical qubits at distance 4 for a 
lattice of radius three. These two labelings are illustrated 
in Figure [Tsl 



C. Discussion 

There are few surprises in our results. For example, 
the best code that we found that maximized the logical 
qubit distance per physical qubit was the compass model 
code in the quadrille tiling, which is already well-known. 
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(a) distance 3 code labelings 




(b) distance 4 code labelings 



FIG. 6: A figure illustrating four labelings of the truncated quadrille tiling (on a radius 1 lattice) that result in a 
quantum code with distance 3/4 (in respectively 6a ' 6b ) , and a number of distance logical qubits proportional to the 

square of the radius of the labeling. 



Furthermore, all of the codes obeyed the upper bounds 
kd G 0{n) and (P G 0{n) — where k is the number of 
logical qubits in the code, d is the distance of the code, 
and n is the number of physical qubits implementing the 
code — that were derived in [30] for codes having spa- 



tially local generators. 

Some of the observed differences between the tilings 
are an artifact of the search space. For example, every 
code found on the hextille tiling could also be imple- 
mented on the deltille tiling, but although we found two 
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FIG. 7: A figure illustrating a labeling of the snub 
quadrille tiling (on a radius 2 lattice) that results with 

distance 4 and a number of distance logical qubits 
proportional to the square of the radius of the labeling. 



kinds of codes for the hextille tiling we found no codes 
for the deltille tilings. This is because our search space 
included no way for the deltille tilings to "knock out" the 
middle qubits in each hexagonal tiling, and furthermore 
the hextille tiling search space included two classes of ver- 
tices which could have independent labelings whereas the 
deltille tiling search space had only one class of vertices. 

Although it is not clear how many of the codes we 
found will have practical applications, the success of this 



search demonstrates the feasibility of using brute-force 
computation to find useful codes within a constrained 
search space. 

IV. CONCLUSION 

In this paper we have presented an algorithm for com- 
puting the optimal quantum subsystem code that can be 
implemented using a given set of measurements. We have 
shown that although this algorithm requires exponential 
time in the worst case, this exponential is a function of 
the code distance, and so the algorithm terminates (rel- 
atively) quickly when the optimal code has low distance. 
Because of this, the algorithm can be used to perform a 
brute-force search through a space of possible measure- 
ments in order to see which give rise to "useful" (high- 
distance) codes. We demonstrated the feasibility of this 
approach by applying the algorithm to search for codes 
implemented on systems with lattice structures corre- 
sponding to nine of the eleven convex vertex-uniform 
tilings, and on all but one of these nine tilings we found 
useful codes. 

This algorithm should prove helpful in two kinds of 
ways in particular. First, it can be applied in an ex- 
ploratory setting to do the tedious work of computing 
the code resulting from a set of measurements so that 
the researcher can experiment with new ideas for choices 
of measurement to see how well they work. Second, it 
can be applied to hone a 'rough' idea for how a code 
might be implemented (such as a particular lattice con- 
figuration) into a concrete idea by scanning through the 
possible choices of the degrees of freedom to see if any re- 
sult in useful codes; of course, cleverness can often come 
up with an answer more quickly than a computationally 
intensive search, but it is good to have the alternative of 
brute-force computation to fall back on when brute-force 
cleverness fails. 
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FIG. 8: A figure illustrating two of the labelings of the snub quadrille tiling (on a radius 2 lattice) that result in a 
quantum code with a single qubit whose distance grows with the radius of the lattice. 




FIG. 9: A figure illustrating two of the labelings of the isosnub quadrille tiling (on a radius 1 lattice) that result in 
a quantum code with a single qubit whose distance grows with the radius of the lattice. 
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(b) distance 4 code labelings 



FIG. 10: A figure illustrating four labelings of the hextille tiling (on a radius 2 lattice) that result in a quantum 
code with distance 3/4 (in respectively 10a' 10b), and a number of distance logical qubits proportional to the square 

of the radius of the labeling. 
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FIG. 11: A figure illustrating two labelings of the truncated hextille tiling (on a radius 1 lattice) that result in a 
quantum code with distance 3 and a number of distance logical qubits proportional to the square of the radius of the 

labeling. 




FIG. 13: A figure illustrating one of the labelings of 
FIG. 12: A figure illustrating two labcHngs of the the hcxadcltillc tifing (on a radius 1 lattice) that results 

truncated hextille tiling (on a radius 1 lattice) that in quantum code whose distance grows with the size of 

result in a quantum code with distance 4 and a number tli^ tiling, 

of distance logical qubits proportional to the square of 
the radius of the labeling. 
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(b) distance 6 code labeling 



FIG. 14: A figure illustrating labelings of the rhombihexadeltille tiling (on a radius 2 lattice) that result in a 
quantum code with distance 5 and four qubits (14a) and a quantum code with distance 6 and two qubits (14b) when 

applied to a radius 3 lattice. 




FIG. 15: A figure illustrating two labclings of the rhombihexadeltille tiling (on a radius 2 lattice) that result in a 
quantum code with distance 4 and sixteen when applied to a radius 3 lattice. 
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